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STRESZCZENIE 

W niniejszej pracy przedstawiono nowy model hydrodynamiczny, ADHYDRO (skrot 

od highly- Anisotropic and strongly-Dissipative HYDRO dynamics), stanowi^cy uogol- 
nienie standardowej hydro dynamiki piynu doskonalego. Umozliwia on efcktywny 
opis zachowania pocz^tkowo silnie anizotropowych ukladow cz^stek produkowanych 
w skrajnie relatywistycznych zderzeniach ci^zkich jonow. Model ten oparty jest 
na przeslankach pochodz^cych z mikroskopowych modeli poczQ,tkowycli stadiow 
ewolucji materii, ktore sugeruj^, ze produkowana materia wykazuje siln^ anizotropi^ 
w przestrzeni p^dow. 

Bazuj^c na ogolnej postaci funkcji rozkladu przeprowadzono analiz^ wlasnosci 
termodynamicznych lokalnie anizotropowych ukladow cz^tek. W wyniku tej anali- 
zy zaproponowano postac uogolnionego rownania stanu umozliwiaj^cego integralny 
opis wielu stadiow ewolucji materii, w ktorym izotropizacja ukladu jest opisywana w 
sposob ci^gly. Nast^pnie wyprowadzono ogoln^ postac rownan opisuj^-cych ewolucji 
silnie anizotropowego plynu. Zaproponowano rowniez postac zrodla entropii zwi^za- 
nego z mechanizmem termalizacji ukladu. W najprostszym przypadku ewolucji jed- 
nowymiarowej w kierunku osi wi^zki oraz przy zalozeniu symetrii ukladu wzgl^dem 
pchni^c Lorentza przedyskutowano szereg wlasnosci modelu. W szczegolnosci, w 
granicy malych anizotropii, wykazano zgodnosc modelu z teori^ Israela-Stewarta. 

Model ADHYDRO po zintegrowaniu z istniej^cym modclem statystyczncj hadroniza- 
cji THERMINATOR zostal uzyty do systematycznej analizy danych eksperymentalnych 
pochodz^cych z eksperymentow prowadzonych na akceleratorze RHIC (Relativistic 
Heavy Ion Collider) w Brookhaven National Laboratory przy energii ^/s^ = 200 
GeV. W ramach modelu uzyskano opis jedno- i dwucz^stkowych obserwabli w zakre- 
sie mi^kkiej fizyki (wartosci p^du poprzecznego czq-stek nicprzckraczaj^ce 3 GeV). W 
szczegolnosci przedyskutowano widma cz^stek w p^dzie poprzecznym, wspolczyn- 
niki przeplywu eliptycznego V2 i przeplywu ukierunkowanego Vi oraz promienie ko- 
relacyjne HBT identycznych pionow. 

Przy zalozeniu istnienia symetrii pchnigc Lorentza w kierunku osi wi^zki w ra- 
mach ADHYDRO wykazano niewielk^ czulosc obserwabli w obszarze centralnych rapid- 
ity na wartosci parametrow charakteryzuj^cych pocz^tkow^ faz^ anizotropow^. Tym 
samym potwierdzono hipotez§ dotycz^c^ wyst^powania uniwersalnego przeplywu 
poprzecznego w zderzeniach ci^zkich jonow bez wzgl^du na anizotropi^ cisnien ten- 
sora energii-pgdu ukladu. 

Ostatccznic, rczygnuj^c z jakichkolwiek symetrii, zastosowano model ADHYDRO w 
najogolniejszej (3+1) wymiarowej postaci do opisu danych poza obszarem central- 
nych rapidity. Ponownie zaobserwowano niewielk^ czulosc obserwabli na szczegoly 
fazy anizotropowej. Bazuj^c na analizie rezultatow modelu uzyskano oszacowanie 
na dlugosc fazy anizotropowej na okolo 1 fm. Wynik ten dopuszcza zatem istnienie 
silnych anizotropii w pocz^tkowych fazach ewolucji materii produkowanej w skrajnie 
relatywistycznych zderzeniach ci^zkich jonow. Jednoczesnie rezultat ten jest zgodny 
z przewidywaniami wielu obecnie stosowanych modeli mikroskopowych. 
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Chapter 1 



Introduction 



1.1 Searching for new states 

of strongly interacting matter 

Relativistic heavy-ion physics is an interdisciphnary research field connecting nu- 
clear physics and elementary particle physics. It emerged in the middle of 1980s, 
when the relativistic beams of heavy ions became available at the Brookhaven Na- 
tional Laboratory and at CERN. The fast development of this field was triggered by 
a broad interest in the experimental and theoretical studies of properties of strongly 
interacting matter in the wide range of thermodynamic parameters such as temper- 
ature, T, and baryon chemical potential, /ib- 

In the framework of Quantum Chromodynamics (QCD), which is the funda- 
mental non-abelian gauge field theory of strong interactions, one expects that the 
strongly interacting matter undergoes deconfinement and chiral phase transitions. 
The experimental verification of those phase transitions is the main aim of the ex- 
perimental heavy-ion programs. 

In the deconfinement phase transition, the hadronic matter changes into a sys- 
tem of interacting quarks and gluons. The initial theoretical concepts, based on the 
phenomenon of the asymptotic freedom (discovered by Gross, Politzer and Wilczek) 
and color screening, indicated that this system was weakly interacting — the asymp- 
totic freedom predicts that the interaction strength between color charged partons 
{quarks and gluons) vanishes at short distances. 

Starting from those ideas, Collins and Perry formulated a hypothesis of exis- 
tence of a new quark state of matter in compressed super-dense star interiors. Also 
at that time, the limiting temperature introduced by Hagedorn was reinterpreted 
by Cabibbo and Parisi as the temperature defining the phase transition between 
hadronic and quark matter. In 1975 they drew the first phase diagram of strongly 
interacting matter. It indicated that at high temperatiirc and/or large baryon den- 
sity the hadronic matter should pass a phase transition to a new state of matter. 
This new phase was termed by Shuryak a quark-gluon plasma (QGP). 

Besides the deconfinement phase transition, at high temperature and/or large 
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baryon density one deals also with the restoration of chiral symmetry, i.e., with the 
chiral phase transition. The order parameter of the chiral phase transition is the 
quark condensate. The standard QCD vacuum is characterized by a non-vanishing 
value of the quark condensate, which is a consequence of spontaneous breaking of 
chiral symmetry at T = /i^ = 0. If the system passes the chiral phase transition, 
the quark condensate drops practically to zero. 

It was quickly realized that the verification of theoretical predictions concern- 
ing properties of strongly interacting matter may be done in laboratory conditions 
by performing collisions of heavy atomic nuclei at ultra-relativistic energies. In 
1980s, various heavy-ion experiments started, collecting data with the required rel- 
ativistic beam energies. First experiments took place in 1986 at the Alternating 
Gradient Synchrotron (AGS) at the Brookhaven National Laboratory (BNL) and 
simultaneously at the Super Proton Synchrotron (SPS) at the European Organiza- 
tion for Nuclear Research (CERN). In 1992 AGS performed ^^^Au+^^^Au collisions 
at the center-of-mass energy per nucleon pair ^snn = 5 GeV. Soon, in 1995 at 
CERN, 208pb^208pb g^g^^g ^ QgY ^g^g collected. The new XXI cen- 

tury started with the construction of the Relativistic Heavy-Ion Collider (RHIC) 
at BNL, which was expected to reach the threshold for the phase transition to the 
QGP phase. Eventually, it has achieved the energy of a/snn = 200 GeV, colliding 
^^''Au beams. In 2010, the largest experimental facility ever built, the Large Hadron 
Collider (LHC) at CERN, designed for colliding ^°^Pb beams, started to collect data 
at the energy of a/snn = 2.76 TeV. 

The experimental results have significantly broadened our general knowledge of 
the properties of strongly interacting matter. First of all, the multi-particle systems 
produced in high-energy collisions of heavy ions turn out to be very much different 
from the systems obtained in more elementary hadron-hadron or electron-positron 
interactions. It is also generally accepted that due to extremely high densities and 
sufficiently large sizes of the systems produced in heavy-ion collisions it is possible 
to ensure the proper conditions for the formation of a quite long-living new state of 
strongly interacting matter. 

Several interesting physical phenomena have been discovered. In particular, 
at RHIC and at the LHC one observes the strong collective behavior of matter 
(quantified by splitting of the slopes of the transverse-momentum spectra of various 
hadron species and by the large elliptic fiow coefficient ^2), strangeness enhancement, 
and strong in-medium jet quenching. These facts indicate that quark and gluon 
degrees of freedom play the main role in the hot and dense medium. However, the 
data suggests also that the produced system is not weakly interacting. The great 
success of relativistic perfect-fluid hydrodynamics in the correct description of the 
data suggests short equilibration times, and consequently large cross sections for 
quark-gluon interactions. This gives hints that the quark-gluon plasma should be 
considered as a strongly interacting fluid rather than as a weakly interacting gas. 

The RHIC experiments suggest the existence of a strongly interacting quark- 
gluon plasma in the high-temperature region of the QCD phase diagram (above the 
critical temperature and for negligible baryon chemical potential, see Fig. II. ip . In 
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Figure 1.1: Phase diagram of the QCD matter. 



this region, the lattice QCD simulations predict a crossover rather than a genuine 
phase transition. In the crossover phase transition, the thermodynamic variables 
change suddenly in a narrow range of temperature, however, all changes are smooth 
and show no discontinuities characteristic for the real phase transitions. The aim 
of the other current (STAR, NA61) and planned (FAIR, NICA) experiments is 
to explore different parts of the QCD phase diagram. The low beam energy scan 
program at STAR and the NA61 experiment search for signatures of the QCD phase 
transition at finite values of the baryon chemical potential. They try to identify the 
position of the critical point that marks the end of a line describing the first order 
deconfinement phase transition at finite values of fiB- The existence of such a line 
is suggested by many QCD inspired theoretical models. 



1.2 Theoretical tools for different stages of 
relativistic heavy-ion collisions 

In order to interpret the multitude of experimental data and formulate conclusions 
about physical properties of matter, we have to compare the experimental results 
with theoretical model predictions. Since, in practice, it is impossible to obtain 
a complete description of such complicated processes within QCD only, we have 
to develop different effective models that are adequate for different stages of the 
evolution. 

Just before the collision, the coherent nuclear wave functions may be viewed (due 
to very high energies) as the gluon sheets described in the theory of color glass con- 
densate (CGC). The passage of the nuclei through each other at the collision proper 
time identified with r = is followed by a subsequent pre-equilibrium stage which 
includes variety of processes leading to the formation of incoherent distributions of 
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partons. 

One argues that the complex configuration of color electric and magnetic fields 
named glasma is produced just after the collision. Its further dynamics is described 
by the classical Yang-Mills theory. Finally, the glasma fields decay by radiating 
gluons. Similar picture is established in the framework of string models where one 
considers collisions of hadrons, followed by the spanning of color strings or tubes. 
In this case, the strings decay producing partons or hadrons. 

The initial dynamics of the collision may be also viewed as reciprocal penetra- 
tion of partons described by their distribution functions, which leads to multiple 
parton scatterings. The subsequent non-equilibrium evolution of the system is then 
described in the framework of kinetic theory with collision terms treated in the 
perturbative limit. This approach is known as the parton cascade model (PCM). 
However, if the system is dense, the particles are off-shell and the use of the quantum 
transport theory becomes necessary (although this is not easy in practice). 

The pre-equilibrium dynamics leads eventually to isotropization of momentum 
distributions of partons, and subsequently to local thermal equilibration of the sys- 
tem. The processes leading to thermalization are not understood well enough at the 
moment, thus the very fast thermalization, expected at RHIC and the LHC, remains 
still an unresolved puzzle. In this Thesis we address this problem and propose how 
to circumvent it. A more detailed discussion of the early thermalization puzzle is 
presented in Chapter [2J 

Multiplicities of particles produced in relativistic heavy-ion collisions are very 
large. Thus, the use of considerations based on statistical physics seems to be well 
justified. If the system reaches the state of local thermal equilibrium, its further evo- 
lution may be described by relativistic hydrodynamics. During hydrodynamic evo- 
lution the local values of thermodynamic variables change gradually. These changes 
are connected with the changes of the underlying phase-space distribution functions. 
If the equilibration rates in the system are not sufficiently fast, the local thermal 
equilibrium is not maintained and the hydrodynamic picture fails. In this case, the 
system's evolution should be described in the framework of the kinetic theory. 

In a series of publications, it has been shown that the space-time evolution 
of matter produced at RHIC may be well described by relativistic hydrodynamics 
of the perfect fiuid P El El H El El H El El IIDI [ID 113 IS] or (if the deviations 
from equilibrium are noticeable but small) by dissipative hydrodynamics with a 
small viscosity to entropy ratio [TH [15]. Relativistic hydrodynamics requires the 
information about the form of equation of state (EOS). The latter can be obtained 
from the lattice QCD simulations (LQCD). The present results obtained at /z^ = 
predict a smooth crossover phase transition from hadronic matter to the quark-gluon 
plasma at the critical temperature Tc ~ 170 MeV. 

The hydrodynamic expansion of the fireball leads to the fast cooling of matter. 
The initially hot and dense quark-gluon plasma passes the phase transition back 
to hadronic matter. The hadronic interactions are less and less effective. At a 
certain stage, the interactions cease and particles freely stream to detectors. This 
process is called the kinetic freeze-out. It is expected that the kinetic freeze-out 
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is preceded by the chemical freeze-out — the stage where hadronic abundances are 
fixed. Nevertheless, in many approaches the identification of the two freeze-outs (the 
single freeze-out model) turns out to be a reasonable approximation. In this case, all 
hadronic observables may be obtained with the help of the Cooper-Frye formula [16]. 
In particular, this approach allows for the successful description of the HBT radii 
consistently with other soft hadronic observables such as the transverse-momentum 
spectra and the elliptic fiow. 

If different hadron species undergo (kinetic) freeze-out at different temperatures, 
the modeling of freeze-out becomes a much more complex problem. In this context, 
there is a growing interest in developing hybrid transport approaches, i.e., the models 
that link hydrodynamics with Monte-Carlo simulations which describe the scattering 
between hadrons in the final state, for example, this is achieved with ultra-relativistic 
quantum molecular dynamics (UrQMD) [T71 118] . 

1.3 Objectives of the Thesis 

In this Thesis we analyze the space-time evolution of matter produced in relativis- 
tic heavy-ion collisions. We propose to describe this evolution in the framework 
of a new, hydrodynamics-like model that inherently includes the possibility of hav- 
ing high anisotropy of pressure at the early stage. We call this new framework 
the ADHYDRO model (highly-Anisotropic and strongly-Dissipative HYDROdynam- 
ics). Using ADHYDRO we try to reproduce the soft hadronic observables measured at 
RHIC. The presence of the anisotropic phase in our description, for about 1 fm/c 
at the beginning of the collision, helps us to circumvent the problem of early ther- 
malization — in order to describe the data we do not have to assume that matter 
reaches local thermal equilibrium within a fraction of a fermi. On the other hand, 
a possible connection of ADHYDRO to the underlying kinetic theory suggests that 
the microscopic collision times are of the order of a few hundredths of a fermi, a 
result suggesting that our picture may be consistent with the notion of a strongly 
interacting quark-gluon plasma. 

The Thesis is organized as follows. In Chapter [2] we address the problem of 
fast thermalization in the very early stages of relativistic heavy-ion collisions. We 
also motivate the development of the ADHYDRO model presented in the Thesis. In 
Chapter [3] we discuss thermodynamic properties of the locally anisotropic systems 
of particles and formulate generalized equation of state. In Chapter H] we derive 
the hydrodynamic equations determining the evolution of matter in ADHYDRO. Sub- 
sequently, in Chapter Owe discuss different aspects of the model in the simple case 
of purely longitudinal and boost-invariant evolution. The Chapter E] contains the 
information about the initial conditions for realistic hydrodynamic simulations. In 
Chapter [7] we include the freeze-out prescription. The results of the model for the 
boost-invariant systems are presented in Chapter [8] and the results for the most 
general non-boost-invariant case are shown in Chapter O We close the Thesis with 
the Summary and Conclusions. 
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Throughout the Thesis we use the natural units where c = h = = I- We 
also use the notation where the metric tensor has the following signature g^" — 
diag(+l, -1,-1,-1). 



The results presented in this Thesis were published in the following papers: 

1. W. Florkowski and R. Ryblewski, 

Dynamics of anisotropic plasma at the early stages of relativistic heavy-ion collisions, 
Acta Phys. Pol. B40 (2009) 2843, arXiv:0901.4653 [nucl-th]. 

2. W. Florkowski and R. Ryblewski, 

Highly- anisotropic and strongly-dissipative hydrodynamics for early stages of rela- 
tivistic heavy-ion collisions, 

Phys. Rev. C83 (2011) 034907, arXiv: 1007.0130 [nucl-th]. 

3. R. Ryblewski and W. Florkowski, 

Non-boost-invariant motion of dissipative and highly anisotropic fluid, 
J. Phys. G38 (2011) 015104, arXiv: 1007.4662 [nucl-th]. 

4. R. Ryblewski and W. Florkowski, 

Highly anisotropic hydrodynamics - discussion of the model assumptions and forms 
of the initial conditions. 

Acta Phys. Pol. B42 (2011) 115, arXiv: 101 1.6213 [nucl-th]. 

5. R. Ryblewski and W. Florkowski, 

Highly- anisotropic and strongly-dissipative hydrodynamics with transverse expan- 
sion, 

Eur. Phys. J. C71 (2011) 1761, arXiv:1103.1260 [nucl-th]. 

6. R. Ryblewski, 

Flow characteristics and strangeness production in the framework of highly- anisotropic 
and strongly-dissipative hydrodynamics, 

to be published in Acta Phys. Pol. B, arXiv:1111.1909 [nucl-th]. 



Chapter 2 

Early thermalization puzzle 



In this Chapter, the early thermalization puzzle in relativistic heavy-ion colhsions 
is discussed. Different theoretical approaches describing early stages of the colli- 
sions are reviewed briefly and different results for the thermalization time-scales 
are presented. Finally, several recent attempts to circumvent the problem of early 
thermalization are introduced. In particular, we provide motivation for the develop- 
ment of our framework of highly-anisotropic and strongly-dissipative hydrodynamics 
(ADHYDRO) which will be presented in greater detail in the next Chapters. 

2.1 Evidence for early thermalization 

The results collected in the heavy-ion experiments at RHIC are interpreted as the ev- 
idence that matter produced in relativistic heavy-ion collisions reaches local thermal 
equilibrium in a very short tim^. The main argument for such early equilibration 
comes from the observation that the space-time evolution of matter is very well 
explained by the relativistic hydrodynamics of perfect fluid [IlEllSllSElinilTllHlini 
\T0\ [11] [T2l IT3] or by viscous hydrodynamics with a small viscosity to entropy ratio 
[m [12] . In particular, a very good description of the soft region of the transverse- 
momentum spectra of hadrons has been achieved in those frameworks. The shape 
of the spectra as well as their azimuthal dependence are interpreted as the evidence 
for the radial and elliptic flow of matter (the latter is described by the elliptic flow 
parameter V2)- The key ingredient of the successful hydrodynamic calculations is 
an early starting time assumed for the hydrodynamic evolution, tq < 1 fm/c. Since 
the initial starting time Tq is identified typically with the thermalization time rth, 
the question arises whether such incredibly fast thermalization of matter can be 
explained by microscopic models of the very early stages of heavy-ion collisions. 

Fast equilibration and perfect-fluidity are naturally explained by the concept 
that the produced matter is a strongly coupled quark-gluon plasma (sQGP) [191 
[2U] . Strong interactions between quarks and gluons might be responsible for the 

""^The physical processes leading to local equilibrium are named by us shortly as "thermalization" 
or "equilibration" processes. 
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fast thermalization rate and small viscosity. Unfortunately, the large value of the 
coupling constant prohibits the use of the standard perturbative methods in this 
case. That is why the methods based on the AdS/CFT correspondence (Anti de 
Sitter/ Conformal Field Theory correspondence) have become now a very attractive 
tool to analyze the early dynamics of heavy-ion collisions. Of course, we have also at 
our disposal LQCD which allows us to study thermodynamic properties of the QCD 
matter pTl [2^2j . In fact, deviations of the lattice results from the Stefan-Boltzmann 
limit in the temperature range T ~ (1 — 3)Tc (where Tc ~ 170 MeV is the critical 
temperature) may be interpreted as the evidence for a strong coupling. 

On the other hand, we should keep in mind that the physics of relativistic heavy- 
ion collisions has been triggered by the idea that at extremely high temperatures, 
T ^ Tc, the asymptotic freedom property of QCD leads to color screening and 
deconfinement [23]. In this case, QGP should behave as a weakly interacting gas 
of quasiparticles (wQGP) treated with perturbative QCD (pQCD) techniques. It 
seems hard to find reasonable explanations for sufficiently fast equilibration of such 
a weakly interacting system. However, several approaches have been proposed in 
order to understand the early time behavior of wQGP. 

The thermalization of wQGP was broadly studied in the framework of Monte- 
Carlo parton cascade models (PCM) where the equilibration of the system was 
understood typically as the effect of multiple hard and semi-hard rescatterings of 
deconfined partons. Originally, PCM included only lowest-order binary (2 — )■ 2) 
collisions [21]. They turned out to be insufficient to provide fast equilibration. 
Moreover, the local equilibration of the system was observed only in the central 
region in coordinate space. Further improvements of PCM models led to inclusion 
of more subtle effects such as: dilated emission from excited partons, the balance 
between radiative emissions and reverse absorption processes, and the interference 
of soft gluons [25l |26l [271 [2H]- They have been proved to significantly affect the 
thermalization rate. Nevertheless, the estimates for equilibration time-scales have 
delivered quite long times, for example, about 1.8 fm/c for Au-|-Au collisions at 
the RHIC top energy. Recently, the importance of multi-particle processes have 
been emphasized. The particle production and absorption processes (2 — )■ 3) [221 
[30] together with the three-particle collisions (3 — t- 3) [31] have turned out to be 
significant in speeding up the equilibration. It has been demonstrated that the 
overall kinetic equilibration may be achieved much earlier but still later than 1 fm/c. 

At the very early stages of heavy-ion collisions, the momentum distribution of 
produced partons is expected to be highly anisotropic. Using theoretical methods 
known from the electron-ion plasma physics in a framework of the quark-gluon trans- 
port theory, one can show that the anisotropies of QGP lead to kinetic instabilities 
with respect to the formation of unstable plasma modes [32l |33]. For QGP, only 
magnetic instabilities known as filamentation or chromo-Weibel instabilities [34j are 
relevant. These strong instabilities drive the system towards the locally isotropic 
state much faster [33 ESI 12Z] than the ordinary perturbative scattering processes 
in the weak coupling regime. This, in turn, indirectly helps to achieve the genuine 
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equilibration 

A quite new and original mechanism of thermalization of QGP has been formu- 
lated during studies of statistical mechanics of the Yang-Mills classical mechanics. 
It has been postulated that equilibration may be understood as the effect of chaotic 
dynamics of the non-Abelian highly excited classical color fields coupled or not to 
the classical colored particles [331 ESI SO]- In this case, the thermalization time 
estimated from the Lyapunov exponent is expected to be in the range 0.4 — 1.0 
fm/c. 

The microscopic models of the early stages of relativistic heavy-ion collisions 
use quite often the concepts of color strings or color-flux tubes. The idea of fast 
thermalization is particularly difficult to reconcile with the results of such models. 
For example, in the color glass condensate (CGC) approach that is often used to 
describe initial stages of relativistic collisions of heavy nuclei [H], the momentum 
distribution functions are far away from the equilibrium ones . 

On the basis of the QCD saturation mechanism incorporated in the CGC model, 
the bottom-up thermalization scenario [Hj was proposed. Since for very large nuclei 
and/or sufficiently high collision energies the saturation scale Qs is much larger 
than the QCD scale Aqcd (at RHIC one expects Qs ~ 1 GeV), the weak coupling 
techniques can be used in the calculations [13]. At the very early times, r ~ Qg^, 
small X gluons with the typical transverse momentum of the saturation scale Qg 
are freed from the incoming nuclei. Subsequently, in the time interval < r < 
a~^^^Qj^, the soft gluons are emitted. After the time r ~ a~^^'^Qj^ the soft gluons 
form a thermal bath and draw the energy from hard gluons. The stop of temperature 
growth identified with full thermalization happens at r ~ a" 

At that time 

the primary hard gluons loose their all energy. Straightforward calculations done in 
this approach for the most central collisions at RHIC give the equilibration time of 
2.6 fm/c or even larger H6] . 

The CGC model itself gives quite interesting picture of the temporal behavior 
of the momentum anisotropies of the matter produced in heavy-ion collisions. For 
the very early proper times, r ^ Q^^, the dominant gluon fields produced in the 
collisions are classical and can be described by the classical Yang-Mills equations. 
In this case, one predicts the energy-momentum tensor with negative longitudinal 
pressure [lH HH]. The pressure asymmetry persists also at later times, r ^ Qj^: 
where the analytic perturbative approaches [19] and the full numerical simulations 
[50] yield the energy-momentum tensor with the longitudinal pressure significantly 
lower than the transverse pressure. This results indicates that the hydrodynamic 
description of the gluon system as a perfect fiuid may start only at later times . 

In the color-fiux-tube models, the production of partons is described by the 
Schwinger tunneling mechanism in strong color fields [511 [52] . The Schwinger mech- 
anism leads to transverse-momentum distributions which have Gaussian shapes. 

^We note that locally isotropic distributions are not necessarily equilibrium distributions. On the 
other hand, the local equilibrium requires that the distributions have the form of Boltzmann dis- 
tributions (or Fermi-Dirac/Bose-Einstein distributions for quantum statistics) which are isotropic 
in the local rest frame. 
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However, if the string tension fluctuates, the Gaussian distributions are changed 
effectively into the exponential distributions [S3] . The resulting forms of the spectra 
resemble the Boltzmann distributions in the transverse-momentum space. Thus, the 
produced system may be viewed as a system where only transverse degrees of free- 
dom are thermalized. This phenomenon may be considered as apparent transverse 
thermalization, since thermal-like spectra are caused by the specific mechanism of 
particle production and not by the parton rescattering processes. 

An intriguing thermalization scenario explaining the equilibration of both the 
transverse and longitudinal degrees of freedom is based on the Hawking- Unruh effect 
[5^ [55] that states that an observer moving with an acceleration a experiences the 
influence of a thermal bath with the effective temperature a/2'K. In the Color Glass 
Condensate theory the typical acceleration in the strong color fields is of the order 
of Qs, which leads to the thermalization time l-n jQs ~ 1 fm/c. 

We conclude our discussion of the early thermalization problem with the state- 
ment that microscopic models of heavy-ion collisions cannot explain naturally the 
thermalization time-scales shorter than 1 fm/c. Moreover, several models suggest 
that thermalization of the transverse degrees of freedom takes place before the com- 
plete three-dimensional thermalization is achievecH. As a consequence, these models 
predict significant momentum anisotropies during early stages of the collisions. It 
means that the standard hydrodynamic picture may be inadequate during very early 
times of heavy-ion collisions. 

These observations have initiated development of effective, hydrodynamics-like 
models which describe the whole evolution of matter and relax the assumption of 
the early thermalization. Typically, this is done by imposing existence of a short 
initial pre-equilibrium phase followed by more standard hydrodynamic evolution. 
These models are discussed below in Section 12.21 

2.2 Relaxing the early thermalization assumption 

Significant differences between the results of various theoretical models of the early 
stages of heavy-ion collisions reflect our lack of precise knowledge concerning the 
processes of particle production. The early stages are, however, of great importance 
as they may have an important impact on further evolution of matter. Based on the 
qualitative hints coming from the microscopic models discussed above in Section [2?H 
we may study different effective scenarios for the early stages. By comparing the 
experimental consequences of such scenarios, we may try to check which hypothesis 
is the most likely to be realized in Nature. 

Of the special interest are the physical scenarios where the assumption of the 
early thermalization of matter is relaxed. In practice, this means that one does 
not start the hydrodynamic description at an early time which is a fraction of a 
fermi but the early stage is described by an effective non-equilibrium model which 



■^Speaking about the complete or full thermalization we have in mind always the approach 
towards local equilibrium. 
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is matched with the hydrodynamic description at a later time. 

In particular, it has been proposed that the early-stage dynamics consists of the 
partonic free-streaming (FS) phase initiated at tq ^ 1 fm/c and followed by a sud- 
den equilibration (SE) process. Typically, the sudden thermahzation is realized by 
applying the Landau matching conditions at a certain transition time rtr > tq. The 
remaining part of the evolution is determined by standard relativistic hydrodynam- 
ics. The first FS+SE approximation was proposed by Kolb, SoUfrank, and Heinz 
[SS] • Recently, it has been reconsidered by several authors [571 EHl EH] • It has been 
demonstrated that the final observables are insensitive to the presence of the early 
free-streaming phase if this phase lasts no longer than 1 fm/c |60] . 

According to the results coming from the string and color-flux-tube models, the 
pre-equilibrium stage may be viewed in a different way, namely, as the transversally 
thermalized parton system. It is conceivable that the initial quark-gluon matter is 
produced already in a state close to the thermodynamic equilibrium in the trans- 
verse direction only, while its longitudinal dynamics is described as free streaming of 
two dimensional clusters [61]. This observation resulted in the development of the 
model that is quite similar to the FS-I-SE framework. In this approach, the initial 
dynamics of the parton system is described by transverse hydrodynamics (TII)[62J. 
Similarly to the FS+SE approach, the system that is initially not completely ther- 
malized undergoes sudden equilibration and its subsequent evolution is described 
by standard hydrodynamics. Interestingly, also in this framework the insensitivity 
of final results with respect to details of the initial transverse-hydrodynamics stage 
has been observed. 

The two approaches discussed above, which have turned out to be successful in 
describing the experimental data, show that the models incorporating high initial 
anisotropy should also account for the fact that matter is eventually fully thermal- 
ized. Therefore, there is a need for models which can, in a comprehensive way, 
describe the gradual transition between early highly-anisotropic stages and later 
fully thermalized stages governed by the perfect-fiuid hydrodynamics. Of course, in 
this context the use of the Landau matching conditions seems to be a very crude 
approximation. 

Certainly, for small deviations from equilibrium the most appropriate description 
should be based on the framework of the second order viscous hydrodynamics [631 
[6^ [65] . However, at the earliest stages of the collision the anisotropics are very large. 
In this regime neither perfect-fiuid hydrodynamics nor viscous hydrodynamics are 
formally applicable. In particular, viscous corrections become compatible with the 
leading terms, which may lead to negative pressure. Of course, these stages can 
be described by models based on the kinetic theory, but the latter has its own well 
known limitations. 

The idea of developing the hydrodynamics-like approach that is suitable for uni- 
form description of several stages of evolution of matter has motivated our devel- 
opment of the ADHYDRO model. This is an effective model which is able to describe 
strongly-dissipative dynamics of initially highly-anisotropic systems. ADHYDRO works 
in three different regimes. At very early stages, it describes the system with ex- 
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tremely high momentum anisotropy. Close to equihbrium, where the anisotropies 
are small, ADHYDRO agrees with the viscous hydrodynamics 0. For very large evolu- 
tion time the system is (locally) isotropic in momentum space and our model gives 
standard perfect-fluid hydrodynamic description. 



*This is strictly valid in the case of purely longitudinal expansion. If the transverse expansion 
is compatible with the longitudinal one, it is possible to switch from the ADHYDRO description to 
the viscous hydrodynamics in the consistent way, as shown in Rcf. .66] . 



Chapter 3 



Locally anisotropic 
systems of particles 



3.1 Locally anisotropic phase-space distribution 
functions 

In this Thesis we consider locally anisotropic systems of particles. We assume that 
such systems are characterized by the anisotropy of their momentum distributions 
analyzed in the local rest frame. The local rest frame at the space-time point x is 
the reference frame in which a small element of the considered system, located at 
this space-time point, is at rest. Locally anisotropic systems of particles may be 
described by the phase-space distribution functions whose dependence on different 
momentum components in the local rest frame is different. 

We restrict our considerations to the case where the anisotropy of the system 
is connected with the difference of pressures in the longitudinal and transverse di- 
rections. These directions are always defined with respect to the collision axis 0. 
The asymmetry in momentum may be realized by introducing two different space- 
time dependent scales, A_l and Ay, that may be interpreted as the transverse and 
longitudinal temperatures, respectively. The explicitly covariant form of such a non- 
equilibrium distribution function is 



/^/(^HW.^). (3.1) 

where p'^ denotes the four-momentum. The four- vector [/^ describes the fluid four- 
velocity 

U>' = j{l,v,,vy,v,), 7 = (l-t;2)-V2. (3.2) 

In comparison with standard equilibrium distributions, we introduce here a new 
four- vector V^. Its appearance is related to the special role played by the beam axis 



^ It is practical to introduce the standard Cartesian coordinate system where the longitudinal 
2-axis is identified with the direction of the beam, while the x and y axes are chosen in such a way 
that the x and z axes span the reaction plane. 
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{z-Sixis) and it is defined in the following way 

V^ = ^M,0,0,1), 7. = (1-^D~'/'- (3.3) 
The four-vectors and satisfy simple normalization conditions 

f/2 = i, V'^ = -l, U-V = 0. (3.4) 

In the local rest frame (LRF) of the fluid element the four-vectors f/^ and have 
simple forms 



^7^^ = (1,0,0,0), \/'^ = (0,0,0,l) 
and the formula (13.11) reduces to 



(3.5) 



f m± 


\P\\\\ 




All ) 



(3.6) 

For massless particles m± = p±, and we may also write 

/ = /(a,ieiii). (3.7) 

Here, for our convenience, we have introduced two dimensionless parameters 

U=P±/X±, eii=P||/A||. (3.8) 

For example, we may consider the Boltzmann distribution that is stretched or 
squeezed in momentum space. 



/o = exp 



exp 



'{p-uy-{p-vy {p-vy 



(3.9) 



where ^ is expressed by the pressure anisotropy parameter x, 

X 



(3.10) 



For isotropic systems x = 1 (^ = 0) and 



/o = fci = exp 



p-U' 



(3.11) 
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3.2 Moments of anisotropic distributions 

The knowledge of the distribution function (13. ip allows us to calculate one-particle 
characteristics of the system under consideration. In particular, the first and second 
moments of the distribution function yield the particle density, energy density, and 
the two pressures. Since at the very beginning of the evolution the considered system 
is expected to consist mainly of gluons, in this Section we restrict our analysis to 
massless particles by setting m = 0. 

Using the standard definitions of the particle number current A^'^ and the energy- 
momentum tensor T^^", 



and 



J (27r)3poP ^ 



(3.12) 
(3.13) 
(3.14) 



and the definition of the entropy current 



where 

$ [/] = Z In ^ - i f 1 + In f 1 + ^ V (3.15) 
fl-o fl-o e V 90/ V 9oJ 

we obtain the following decompositions (the method leading to Eqs. fl3.16l) - fl3.18l) is 
explained in Ref. p7]): 

= nU", (3.16) 
T'^^ = {e + Px) Uf'U" - Pi_ g"" - (P^ - P||)^^^^ (3.17) 

and 

S'' = aU''. (3.18) 

Here Qq is the degeneracy factor related to internal quantum numbers of particles 
that form the fluid and = Ep= ^m? + p\ + p| • For simplicity, in Eq. (13.151) we 

have introduced e = +1 for bosons and e = — 1 for fermions. In the limit e — we 
obtain the classical Boltzmann definition. 

Equations fl3.16l) - fl3.18l) contain non- equilibrium parameters: particle density n, 
entropy density a, energy density e, transverse pressure Pi, and longitudinal pres- 
sure P|| E|. 

^Below, the parameters such as n, a, e, Pi and P|| wih be cahed and treated as thermodynamics- 
hke or thermodynamic parameters, although, strictly speaking, these quantities do not describe 
the equilibrium state. 
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In the isotropic case, the longitudinal and transverse pressures are equal, P± = 
P\\ = Peq, and we recover the standard form of the energy-momentum tensor valid 
for the perfect fluid. In general, when the system is out of equilibrium, the last term 
in f l3.17p is different from zero. Then, in the local rest frame of the fluid element we 
findEn 

/ £ \ 
P± 
P± 
\ P|| / 



(3.19) 



Hence, as expected, the structure (13.171) allows for different pressures in the lon- 
gitudinal and transverse directions. For massless particles the condition that the 
energy-momentum tensor is traceless leads us to the following expression 



2P, +P 



(3.20) 



From Eq. (13.201) we can infer that thermodynamics-like parameters describing the 
anisotropic system are functions of two independent parameters. For example, the 
energy density is a function of two pressures, e = e{P±,P\\), or equivalently of two 
temperatures, e = £:(A_l,A||). This feature distinguishes anisotropic systems from 
equilibrated matter where the energy density for baryon-free matter is a function 
of isotropic pressure only, e = e{Peq). For similar out-of-equilibrium systems with 
non-vanishing baryon number, three parameters are needed to specify the local state 
of the system. 



3.3 Pressure anisotropy 



The internal consistency of the approach based on the distribution function (13.11) 
allows us to derive the following concise structure of the thermodjTiamic parameters 
describing the anisotropic system of particles (for more details see |67] ) 



e{x, n) 
P±{x,n) 
P\\ (x, n) 



4/3 



- ) Rix), 

9/ 

9 J 
n\ 
9 J 



4/3 



R{x) 



+ xR'(x) 



2xR'(x) 



where the function R{x) is defined by the integral 



27r2 



(3.21) 
(3.22) 
(3.23) 



(3.24) 



^The Lorentz transformation leading to LRF is defined explicitly in 
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and g is a constant defined by the expression 

S = l'-^^f«^M (3.26) 

The range of the integration over and is always between and infinity. The 
function R'{x) in Eqs. fl3.22p and (13.231) is the derivative of R{x) with respect to x. 

The pressure relaxation function R{x) describes bulk properties of the system 
characterized by the anisotropy x. The particular form of R{x) follows directly 
from the microscopic calculation if the distribution function (13. ip is known. At the 
more phenomenological level of description, a special form of the function R{x) can 
be assumed and treated as an ansatz describing the anisotropic system [69]. The 
anisotropy parameter x and the non-equilibrium entropy density a are the most 
convenient independent variables used in our further analysis. 

Equations (I3.2ip - (l3.23p impose general physical constraints on the form of R{x). 
For X <C 1 (A_L <^ A||) the function R{x) behaves like x~^^^. In this limit P± = 
and e = P\\. Similarly, for a; 1 (A_l ^ Ay) the function R{x) behaves like x^^^, 
implying that Py = and e = 2P±. This behavior is expected if we interpret the 
parameters A_l and Ay as the transverse and longitudinal temperatures, respectively. 
Since color fields are neglected in this approach both pressures are expected to be 
positive. Pi > and Py > 0. Furthermore, close to equilibrium the longitudinal 
and transverse pressures are expected to equilibrate, Pj_ = Py. Thus, according to 
Eqs. (I3.22P and (13.231) we obtain another key restriction for R{x), namely, 

P'(l) = 0. (3.26) 

Based on the RHIC experimental data, we assume that the system eventually 
achieves the perfect fluid limit, hence we may write that lim a; = 1. 

T— >00 

The particle and entropy densities introduced in Eqs. (13.160 and (13.180 are de- 
fined by the following integrals 

oo oo 

n = ^jdUU j di\\f{U,i\^, (3.27) 



oo oo 

^ = / ^^^^^ / ^^11 ^ (^^'^ll)] • (3-28) 



Comparison of Eqs. (13.271) and (I3.28P indicates that the particle density and the 
entropy density are proportional to each other, with the proportionality constant 
depending on the specific choice of the parton distribution function /. 

3.4 Special forms of anisotropic distribution func- 
tions 

So far we have not considered any explicit forms of the distribution functions. In this 
Section, we turn to the discussion of the anisotropic distribution functions satisfying 
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the ansatz (13. ip and having the following form 



f = 9o 



exp(^ei + eS) -e]"'. (3.29) 



Equation fl3.29p may be regarded as a simple generalization of the equilibrium dis- 
tributions in the case of vanishing baryon number. In a more general case, we should 
take into account also the net baryon number density ns and introduce the baryon 
chemical potential hb in fl3.29p . However, if we focus on the central rapidity region 
in the ultra-relativistic limit (which is the case when analyzing the RHIC and LHC 
experimental data) we may neglect the baryon chemical potential fiB- Hence, the 
thermodynamic parameters remain the functions of entropy density and anisotropy. 

The function /(e = 0) may be considered as the Boltzmann equilibrium distri- 
bution that has been stretched (or squeezed) in the longitudinal direction. For com- 
pleteness, it is also interesting to consider the stretched (or squeezed) Bose-Einstein 
distribution /(e = +1) = /be or Fermi-Dirac distribution, /(e = —1) = fpo- We 
note that in the local equilibrium limit, Eq. (13.291) with e = is reduced to the 
standard Jiittner formula where A_l = Ay = T. 

With the definition (K^ . we recover the structure flX^ - flX^ P]. In the 
case of the classical Boltzmann statistics, the relaxation function (13.240 is given by 
the form El 



and the constant fl3.25p is simply 



X arctan vx 
1 + - 



(3.30) 



g = ^. (3.31) 

Equations (13.270 and (I3.28p relate the entropy and particle densities in the following 

way 

a = An. (3.32) 

In the similar way we may consider the Bose-Einstein and Fermi-Dirac distributions, 
see summary in Table 13.11 . 

The results presented in Table [XT] may be used in Eqs. fl3.2ip -( l3.23p . In this 
way we find that the numerical factor e/ebe (giving the ratio of the classical energy 
density to the Bose-Einstein energy density for the fixed a and x) is very close to 
unity, e/ebe ~ 1.027 . It means that the equations defining the energy density and 
pressures in terms of a and x for the classical and quantum systems are practically 
the same. Therefore, in our further studies we use the classical statistics given by 
Eq. (13.291) . Moreover, in order to account for the internal degrees of freedom of 
gluons we set qq = IQ (spin x color). 

We stress that physical restrictions concerning the pressure relaxation function, 
discussed in Section 13.31 exclude from the consideration certain forms of distribu- 
tion functions that formally satisfy the ansatz fl3.ip . For example, the anisotropic 



''Note that for a: < 1 the function (arctan — l)/\/x — 1 should be replaced by 
(arctanhA/1 — x) / a/1 ~ x 
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classical particles bosons fermions 

9 f C(3)^ C(3) (1-2-2)^ 

a{n) a = An a = a = 



Table 3.1: Thermodynamic properties of the anisotropic systems for different statistics. 
The rows contain the following information: R{x) — pressure relaxation function ()3.24p . 
g — constant defined by Eq. (j3.25p . cT(n) — relation between entropy density and particle 
density calculated from Eqs. (|3.27p and ()3.28p . 

distribution functions of the factorized form /2 = goexp {—$,±) exp do not 

satisfy the condition fl3.26p . 

In the upper part of Fig. I3.1l we present the rescaled pressure relaxation function 
r{x) = 7c'^R{x) I (3(70) while in the lower part of Fig. 13. II we have plotted its derivative 
r\x). We can see that according to (13.261) r'[x) vanishes when Pj_ = Py. In Fig. 13.21 
we show the ratios: P\\IPl (solid red line), P\\Ie (green dotted line), and P^je (blue 
dashed line) for the case of R{x) given by Eq. (13.301) . The considered ratios are 
functions of x only. In agreement with previously given remarks we see that e = P\^ 
for a; = and e = 2P± for x — ?■ oo. The equilibration of pressures is reached for 
a; = 1. To a good approximation P\\/P± ~ x~^^^, thus, the anisotropy parameter x 
may be treated as the direct measure of the pressure anisotropy. 



3.5 Generalized equation of state 

The equations of perfect-fiuid hydrodynamics form a closed system of equations only 
if they are supplemented with the equation of state that specifies thermodynamic 
quantities as functions of one parameter [§. Similarly, our framework requires that all 
thermodynamic quantities should be expressed by two arbitrarily chosen parameters, 
e.g., by the non-equilibrium entropy density a and the anisotropy parameter x. Such 
relations play a role of the generalized equation of state and allow us to close the 
system of dynamic equations (the latter will be formulated in Chapter H]). 

If we consider a system of partons described by the anisotropic phase-space 
distribution function given by Eq. (13.291) and obtained by squeezing (or stretching) 
the classical Boltzmann distribution, we have to use Eqs. (I3.21I) - (I3.23I) together 
with Eqs. (I3.30p - (l3.32p . Equations (I3.2ip - (l3.23p describe the non-equilibrium state 
defined by the values of the non-equilibrium entropy density a and the anisotropy 
parameter x. Equations (I3.21I) - (I3.23I) approach naturally the equilibrium limit if 



^If the baryon-free matter is considered. 
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Figure 3.1: a) The function r(x) (solid red line) shown together with the asymptotics 
r(x) ~ {'k/A)x^/^ (blue dashed line) and r(x) ~ [1/2) x^^^^ (green dotted line), b) 
The derivative r'(x) of the pressure relaxation function (solid red line). 



X — )■ 1, since r(l) = 1 and r'(l) = 0. In the equilibrium the following relations are 
vahd: eid = 3^oTV7^^ Pid = c/oTVvr^ and did = ^o^Vvr^ [70j. Hence, Eqs. 
(13.231) may be written equivalently in the form, 



e{x, cr) 
Pj_(x,cr) 
P||(a;,a) 



£:id(a)r(x), 

-Pid(o") [r{x) + 3xr'(x)] , 
Pid{,(^) [r{x) — 6xr'(x)] , 



(3.33) 



where ei^{a) and Pid(o") are equilibrium expressions for the energy density and pres- 
sure. 

Equation of state fl3.33p is very likely to be realized at the very early stages 
of the collisions, when the produced system is highly anisotropic and described 
by the distribution functions discussed in this Chapter. On the other hand, as 
the system expands, the interactions between its constituents thermalize it, so the 
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100 



Figure 3.2: The ratios P\\/P± (solid red line), P\i/e (green dotted line), and P±/e 
(blue dashed line) shown as functions of the variable x, the results for the pressure 
relaxation function fl3.30p . We also show approximation of P\\/P± by the function 
3;-3/4 (^i^iaci^; dot-dashed line). 



system's thermodynamic variables approach the realistic EOS and x tends to unity. 
The realistic EOS for the strongly interacting matter, encoded in the temperature 
dependent sound velocity, is shown in Fig. 13.31 At low temperatures it is consistent 
with the hadron-gas model with all experimentally identified resonances. On the 
other hand, at high temperatures it coincides with the LQCD results. 

The two limiting cases discussed above suggest that one may use the generalized 
equation of state of the form 



-qgp 



»r(x). 



(3.34) 



P±{x,(t) 
P\\{x,a) 



^qgp(cr) [r{x) + 3xr'{x)] , 
-Pqgp(o') [r{x) - 6xr'{x)] . 



Here, the equation of state of the ideal relativistic gas, defined by the functions £id(c") 
and Pid(cr) has been replaced by the formulas characterizing the realistic equation 
of state for vanishing baryon chemical potential, i.e., by the functions eqgp(o") and 
-Pqgp(cr) corresponding to the sound velocity function depicted in Fig. 13.31 

We have to stress that there is no microscopic explanation for the form of 
Eqs. (13.341) . The main reason for this is that we cannot provide simple quasi-particle 
picture for the realistic QCD equation of state. Despite this Eqs. (13.341) have many 
attractive features. Firstly, as discussed above, they reproduce the two attractive 
limits and they interpolate between them in the continuous way. Moreover there is 
no additional assumptions needed (such as, e.g., the Landau matching conditions) in 
order to describe the change from one thermodynamic regime to another one. Ad- 
ditionally, Eqs. fl3.34l) naturally include the phase transition from the QGP to the 
hadron gas. In this way Eqs. (I3.34p are able to describe different stages of evolution 
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Figure 3.3: Temperature dependence of the square of the sound velocity for strongly 
interacting baryon-free matter obtained in Ref . [7l] . The results of the lattice simulations 
of QCD are represented by the solid line, and the results obtained with the hadron- 
gas model are represented by the dashed line. The thick solid line describes a simple 
interpolation between these two calculations in the temperature region close to Tc = 170 
MeV. 



of matter in heavy-ion collisions in a uniform way. 



Chapter 4 



Highly-anisotropic hydrodynamics 
with strong dissipation 



In this Chapter we introduce the dynamic equations that determine the evolution of 
a highly-anisotropic fluid. Our approach is based on the energy-momentum conser- 
vation law and on the postulate that specifies the entropy production in the system. 
The resulting equations define the ADHYDRO model. 

4.1 Space-time dynamics 

of highly-anisotropic systems 

In Chapter [3] we argued that the system formed at the very early stages of relativistic 
heavy-ion collision may be described by the anisotropic phase-space distribution 
function (13. ip . We assume now that the space-time dynamics of such a system is 
governed by the following equations, 



Equation fl4.ip expresses the energy and momentum conservation law. It has an 
analogous form to the equation that defines dynamics of the perfect fluid. On the 
other hand, Eq. (14. 2 p describes the entropy production determined by the entropy 
source S. In the standard perfect-fluid approach, the right-hand-side of Eq. (14.21) 
vanishes and Eq. (14. 2 p becomes the condition of adiabaticity of the flow. We note 
that in the perfect-fluid hydrodynamics Eq. (14.21) follows directly from Eq. (14. ip and 
the appropriate form of the energy-momentum tensor. 

In our approach, Eq. (14. 2 p should be treated as an ansatz. In Chapter [3] we 
showed that the energy-momentum tensor T^^" and the entropy flux appearing in 
Eqs. (14. 1 1) - (14. 2 p have the structure defined by Eqs. (I3.17P and (13.181) . In this case, 
the specific form of the entropy source S allows us to close the system of equations 
and determines the dynamics of the anisotropic fluid 0. 

-"^This issue will be discussed in more detail later in Section [131 



0, 



(4.1) 
(4.2) 
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The projections of Eq. fl4.1l) on the four-vectors Ui, and V^, yield 

W'd^e = -{e + P^)d^U^+{P^-P\\)U,V^d^V'', (4.3) 
V^'d^Pw = -{P\\-P±)d,y^ + {e + P±)VM^d^U\ (4.4) 

Equation (14. 3 p expresses the energy conservation, whereas Eq. (14. 4p describes the 
conservation of the longitudinal momentum. 

4.2 General three-dimensional expansion 

In the general case, where matter expands in the longitudinal and transverse direc- 
tions without any symmetry constraints, we may use the following parametrization 
of the four- velocity of the fluid If^ and the four- vector 

= {uq cosh {},Ux,Uy,Uo smb. (4.5) 
V^" = (sinh ^9, 0,0, cosh ^9). (4.6) 

In Eqs. (14. 5p and (14. 6 p we have introduced the two components of the transverse 
four-velocity of the fluid, and Uy, and the longitudinal fluid rapidity Using the 
normalization condition (13. 4p we also find 



Mo = \/l + ^^i, u± = Jul + ul- (4.7) 



If the parametrizations (14.50 and (14. 6 p are used, the differential operators appearing 
in (14. 3 P and (14.41) take the following form 







(4.8) 




= V_L ■ U_L + LiUo + UqL2^, 


(4.9) 




= L2, 


(4.10) 




= L.'d, 


(4.11) 






(4.12) 




= -Uq (u_l ■ V_L + UqLi) 


(4.13) 



where the two linear differential operators Li and L2 are defined by the expressions 

d 

Li = coshy^^-sinhy^, (4.14) 

r 

d 

L2 = -sinhra^ + coshr^. (4.15) 

r 

In Eqs. (I4.14p - (l4.15p we have introduced the space-time rapidity 77 and the proper 
time r, 

1 , t + z , , 



r 



Vt^ - z^. (4.17) 
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The difference between the longitudinal ffuid rapidity and the space-time rapidity 
T] is defined as 

¥ = 7]-^. (4.18) 

We use the boldface notation for two-dimensional vectors in the transverse plane: 
u± = {ux,Uy) and V± = {dx,dy). 

Besides Eqs. f l4.3p and (14. 4p . the two additional equations describing transverse 
dynamics of the system are needed. They can be chosen, for example, as the two 
linear combinations: Uid^T^^ + U2df,T^^ = and U2d^T^^ - Uid^T^""^ = 0. In this 
case we find 



VU I : 



e + Pi 



' ^^^^ + VP^ + (P^ - Pn)U,V^d^V' 
u I 



(4.19) 



Here = U^d^ is the total time derivative and 6 = d^U^ is the volume expansion 
rate. Using this notation, the entropy production equation, see Eq. (14. 2p . can be 
put in the compact form 

Va + ae = T.. (4.21) 

Equation together with Eqs. (gJD, (|13D, fICTD and fICTD form a set 

of hydrodynamic equations describing dynamics of the anisotropic fluid. We refer 
shortly to this system of equations as to the ADHYDRO equations or as to the ADHYDRO 
model or framework. 

If the generalized equation of state e{a,x), see Sect. 13.51 and the entropy source 
S((j, x), see Sect. 14. 3[ are specified, the ADHYDRO equations form a closed system 
of five equations for five unknown functions: two components of the fluid velocity 
u^; and Uy, the longitudinal rapidity of the fluid the non-equilibrium entropy 
density a and the anisotropy parameter x. These functions depend on transverse 
coordinates x, y, the space-time rapidity r], and the proper time r. One has to solve 
the ADHYDRO equations numerically for a given initial condition specified at a certain 
initial time tq (to be discussed in Chapter [6]) . 



4.3 Entropy source 

The entropy source S describes the entropy growth due to equilibration of pressures 
in the system. As we have already mentioned above, Eqs. (14. ip and (14.21) may be 
solved only if the function S((T, x) is specified. The functional form S(cr, x) must be 
delivered as the external input for the anisotropic hydrodynamics. 

The simplest ansatz for S((T, x) satisfying general physical assumptions may be 
proposed in the following form [68l [72|, 173] 



(4.22) 
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where Teq is a time-scale parameter introduced to control the rate of the processes 
leading to equilibration of the system. The expression on the right-hand-side of 
Eq. fl4.22p has several appealing features. First of all it is positive, as expected 
on the grounds of the second law of thermodynamics. Moreover, it has a correct 
dimension, since S is proportional to the entropy density a. In this way Eq. (14.221) 
does not destroy the scale invariance of the perfect fluid hydrodynamics, which 
allows for multiplication of a in the evolution equations by an arbitrary constant. 
In the natural way, the entropy source S defined by (14.221) vanishes in equilibrium, 
where a; = 1, thus Eqs. (I4.ip - (l4.2p can be reduced to the equations of perfect 
fluid hydrodynamics in the limit a; — ?■ 1. The expression (I4.22p is also symmetric 
with respect to the interchange of A_l with Ay (consequently S does not change if 
X ^ 1/x). Finally, for small deviations from equilibrium, where |a; — 1| -C 1, we find 

(x - 

So(x) ^ ^—-^a. (4.23) 

The quadratic dependence of the entropy source displayed in (I4.23P is characteristic 
for the 2nd order viscous hydrodynamics (to be discussed in Chapter 15.31) . 



4.4 Anisotropy evolution 

If both the generalized equation of state and the entropy production term are de- 
fined, one can derive a compact expression for the time evolution of the anisotropy 
parameter x. Of course, this expression should be considered together with other 
dynamic equations. Nevertheless, its form turns out to be useful in the analysis of 
the time evolution of the system. 

In the case described by Eqs. (I3.34p we use Eq. (14.31) and find 

Vx = ^^^(3U,V^d^V-9)-(l + ^]^, (4.24) 



^qgp V ^qgP / '^eq 



where 

H{X) = ^-Teq. (4.25) 

rix) a 



Similarly, in the case described by Eqs. (I3.33P we find 



Vx = xiSU^V^d.V -9) - ^H{x). (4.26) 

3r„, 



'cq 



Equation (14.261) may be treated as the special case of (14.241) where we 



qgp- 



Chapter 5 



Purely longitudinal 
boost-invariant motion 



In this Chapter we study a simple case of purely longitudinal, boost-invariant ex- 
pansion. This situation is easy to investigate and we can obtain useful information 
about the behavior of anisotropic systems by applying simple analytic methods. 
This will help us to understand the behavior of more complicated systems which 
either expand transversally or do not exhibit boost invariance. 

5.1 Implementation of boost invariance 

The measurements done by the BRAHMS Collaboration [71] at RHIC showed that 
the rapidity distributions of charged particles are almost constant in the midrapid- 
ity region (|y| < 1). This observation suggests that all measured quantities, to a 
good approximation, are boost-invariant in this region. The boost-invariance implies 
that all thermodynamic variables are functions of the proper time r and transverse 
coordinates only, while the longitudinal flow has the form 



Since t = r cosh 77 and z = rsinh?], Eq. fl5.ll) is equivalent to the simple condition 



i.e., Y = 0. With this assumption at hand, it is sufficient to find the solutions of 
our model for z = {t = t), since solutions for other values of z may be obtained 
by simple Lorentz transformations. 

To illustrate the main features of our model with simple examples, we neglect 
now the transverse expansion of the fiuid. This is a reasonable approximation 
for the early stages of the evolution, where the longitudinal expansion dominates, 
Ux = Uy = {uq = 1). We also assume now that the system may be treated as ho- 
mogeneous in the transverse plane, hence, we neglect dependence of thermodynamic 



Vz = tanh?? = -. 



(5.1) 



(5.2) 
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quantities on the transverse variables x and y. These assumptions are characteristic 
for the Bjorken model [75], where the behavior of the perfect fluid has been consid- 
ered. In our approach the perfect-fluid dynamics is replaced by the hydrodynamics 
of the highly-anisotropic fluid. 

In the considered one-dimensional and boost-invariant case Eq. (14. 4p is auto- 
matically fulfilled, since V = dr, 6 = 1/t, and rV'^d^V'^ = W^. After neglecting 
transverse expansion, also Eqs. (14.191) and (14.201) are satisfied automatically. There- 
fore, we are left with Eqs. (14. 3 1) and (14.211) which are reduced to the expressions 

^ - (5.3) 

dr r 

rfa 1 S , ^, 

-- + - = -. 5.4 
aar r cr 

Equations (15.31) and (15.41) are two equations for two unknown functions, (j(r) and 
x(r), which may be solved numerically. 



5.2 Asymptotic limits 

As we have discussed above in Sect. 14.41 Eqs. (15. 3p and (15. 4p may be used to de- 
rive an equation describing the evolution of pressure anisotropy x. For the purely 
longitudinal boost-invariant motion, Eq. (I4.26P simplifies to an ordinary differential 
equation for x, 

dx 2x AH{x) 
dr T 3req 
which with the help of Eq. (I4.25P may be also rewritten as 

-'(^) - I) = • (5.6) 

First of all, it is useful to consider the isentropic flow where S(cr, x) = 0. In this 
case, from Eq. (15. 4p we obtain the Bjorken solution 

a= , (5.7) 

r 

where the parameters ctq and Tq are arbitrary constants. For S(cr, x) = the right- 
hand-side of Eq. (15.61) vanishes. This implies that either x = 1 oi x = xqt'^/tq 
(where xq is another integration constant). 

The case x = 1 corresponds to the standard, perfect-fluid hydrodynamics. In- 
deed, for an ideal relativistic gas in equilibrium we have P\\ = P± = Peq and Eq. (15.31) 
implies e(r) ~ t~^^^. In the second case we obtain a simple longitudinal free- 
streaming solution (see our detailed discussion in Section 15.51) . This situation is 
equivalent to the case where the equilibration time is set to infinity, Teq = oo. In 
addition we can consider also the case where Py = 0. In this situation Eq. (15.31) 
describes the evolution of the energy density valid for transverse hydrodynamics 
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5.3 Consistency with Israel-Stewart theory 

For a near-equilibrium system, the distribution function can be decomposed into an 
equilibrium part plus a small correction part. The methodology of including such 
corrections in the consistent way leads to the framework of relativistic viscous hy- 
drodynamics of the first order (non-causal Navier-Stokes theory) or the second order 
(causal Israel- Stewart formalism [531 Ell When deviations from equilibrium are 
small, which in our approach may be quantified by the condition |,^| ^ 1, and the 
motion is purely longitudinal, one can show that our model is consistent with the 
Israel-Stewart theory. 

As before, we consider again the baryon-free systems. The 2nd order dissipative 
fluid is characterized by the energy-momentum tensor 

where the pressure tensor P^^ is defined by the formula 

p^^- = -(p^q + n)A'^" + TT^". (5.9) 

Here A'^" = g^'^ — U^V^ is the tensor projecting on the space perpendicular to the 
four-velocity f/^, Peq is the isotropic pressure, 11 is the viscous hulk pressure, and 
T^^^ is the viscous shear tensor. 

The energy-momentum tensor (13.171) . which forms the basis of our approach, 
may be rewritten in the form (15.81) if we make the following identification 

pi^u ^ .p^A^'^ - (Pi - P\\)V^V''. (5.10) 

Let us ignore now the bulk viscosity, since for massless particles or ultra-relativistic 
particles the bulk viscosity vanishes, 11 = 0. If there is no transverse expansion, 
we may compare Eqs. (15. 9p and (15.101) in the Landau local rest frame, and make 
the identification of the longitudinal and transverse pressures: P±_ = Peq + 7f/2, 
P|l = Peq — 7f, where Peq = ^eq/S = e/3 and vf = vr^^ is the 33 component of the fluid 
shear tensor. Then, the expansion for small ^ gives flEl [77] 

eeq R{x) 45^ 
By differentiating Eq. (I5.1ip with respect to the proper time we obtain 

4:5 f 1 dW W de\ 

Y [id^ ~ l^d^ J^d^' ^^'^^^ 

Incorporating the leading order terms in ^ of Eqs. (15.31) and (15. 5p in Eq. (I5.12p one 
finds 

dif 4:W 16 e 15W 

— H = . (5.13) 

dr 3t 45 t 4req ^ ^ 
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We can match our model with Israel-Stewart theory by making the following iden- 
tifications 

1 4 5r]^ 



(5.14) 



where x^- is the shear relaxation time, rj.^ is the shear viscosity, T is the temperature 
in the isotropic system, and cxeq is the equilibrium entropy density. In this way, 
Eq. fl5.13p takes the form of the 2nd order viscous hydrodynamic equation for vf 

dW Aw AtIt^ 7f 

^ = -37+3^--. (5.15) 
The entropy production in the Israel-Stewart theory is given by the formula 



d.S^ = (5.16) 



and for the longitudinal expansion considered in this Section, Tiaj3T^"'^ = 37r^/2. 
Substituting Eqs. (15. lip and (15.141) into Eq. (I5.16P we obtain 



?2 



d.S'^ = aeq-^, (5.17) 

which is consistent with our main ansatz for the entropy production if deviations 
from equilibrium are small. 

One should stress that deviations from equilibrium in the Israel- Stewart theory 
have to be small, so that the dissipative fiuxes are small compared to equilibrium 
variables, y/n^^nt^'^ <^ Peq- If this condition is not fulfilled, the use of the Israel- 
Stewart framework cannot be formally justified. 



5.4 Deriving entropy production 
from kinetic theory 

Our formulation of the highly-anisotropic hydrodynamics in one-dimensional case 
was followed independently by Martinez and Strickland [771 [75] , who derived equa- 
tions analogous to Eqs. (I5.3l) - (l5.4p starting from the relativistic Boltzmann equation 

p>^dJ{x,p) = C[f{x,p)]. (5.18) 

Here C[f{x,p)] is the collision kernel which accounts for all microscopic interaction 
processes among particles. In the Martinez-Strickland model [771 178] , the evolution 
of matter is described by two coupled ordinary differential equations: the first one 
for the momentum anisotropy ^, and the second one for the typical hard momentum 
scale Phard (average momentum in the parton distribution function). These equations 
are obtained by taking the 0th and 1st moments of the Boltzmann equation (I5.18p . 
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In general, the entropy production may be calculated directly from the definition 
(13 .141) with the help of Eq. fl5.18p . For classical particles, it is given by the expression 



(27r)V 



-C[/(x,p)](ln/(x,p) + l). 



(5.19) 



Martinez and Strickland used the relaxation time approximation for the collision 
kernel. This allows for a direct calculation of the right-hand-side of Eq. f l5.19p and 
gives the entropy source Sm5 (analogous to our S), 



+ 



RusiO ' 3(1 + 



(5.20) 



MS 
eq ■ 



In Eq. f l5.20p . the parameter T denotes the inverse of the relaxation time, F = l/r^ 
The function Hms{^{x)) appearing in Eq. fl5.20l) is an analog of our function H{x 

3(e + 1) (v/e+T<^s(0 - l) ^Ms(0 



^^^^^^ 2Rusio + m + m'usiO 

and Rms{C{^)) is an analog of our R{x), 



(5.21) 



1 arctan a/? 
+ 



1 + e 



vr 

36-0 



[l+^)-y^R{l+^). 



(5.22) 



For small anisotropies, one can check that the Martinez-Strickland approach coin- 
cides with our model if 

' eq ^ ' cq 



5.5 Longitudinal free-streaming 

In the case of anisotropic distribution functions (see Section 13. ip it is interesting 
to consider the pure longitudinal free transport (alias free- streaming) of particles, 
which is the regime in which particles do not interact completely with each other. 
In this case, the distribution functions should satisfy the coUisionless Boltzmann 
kinetic equation 

p^djix,p)=0. (5.24) 

We assume that the effects of free-streaming do not change the shape of the anisotropic 
distribution function fl3.ll) and the time evolution of the parameters A^, Ay, and 
is determined by the conservation laws. 

We parametrize the four-momentum of a particle in the standard way 



pf" = (Ep,p_L,P||) = (m_LCOshy,p_L,m_Lsinhy) , 



(5.25) 
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where we have introduced the transverse mass m_|_ 
a particle 



y 



2 Ep- 



+ and the rapidity of 



(5.26) 



For simphcity, in the following analysis we assume that particles are massless, m = 0. 
For the pure longitudinal boost-invariant expansion we use Eqs. (13.11) . (14. 5p . (14. 6p . 
and (15.251) . and we rewrite Eq. (15.241) in the form 



d sinh(y - rj) d 

cosh(y -VJtT^ ^ 

or T or] 



f{w,v)=0, 



(5.27) 



where we have introduced the variable w = p±/X±{r) and v = p± sinh(y — r])/ Ay (r). 
After simple algebra we obtain 



df dX± ^ sinh(y — r]) df 
dw dr 



+ 



dv 



dr 



All 



0. 



(5.2^ 



This equation is satisfied by any form of the function / provided A_l = A^ and 
A|| = Ajjro/r, where Tq, A*} and Ay are constants. In the considered case the anisotropy 
parameter x has the form 

' AO 



X 



Ajjro 



(5.29) 



where xq = specifies the initial anisotropy of the system. Thus, we see 

that our approach with Teq — ?• oo includes the boost-invariant free-streaming as a 
special case. 



5.6 Results for purely longitudinal boost-invariant 
motion 

We solve Eqs. (15.31) and (15. 4p numerically with the entropy source defined by Eq. (I4.22p 
with the initial conditions x = xq and a = ao specified at the initial proper time 
T = Tq. For comparison, we also use the entropy source defined by Eq. (15.201) . The 
time-scale parameter Teq is set equal to 0, 0.25 fm, 0.5 fm, 1 fm, and oco. The values 
of the relaxation time and the parameter F are connected by Eq. (15.231) . In all 
the numerical calculations, the initial conditions are specified at the initial proper 
time To = 0.25 fm. We consider three cases: 

i) Xo = 100, where the transverse pressure dominates over the longitudinal pres- 
sure during the very early stages of the evolution (as suggested by the micro- 
scopic models), in this case the momentum shape is oblate, i.e., it is stretched 
in the transverse direction and squeezed in the longitudinal direction, 

^The value Teq = means that in the numerical calculations we take Toq = 10^'^ fm, similarly, 
the value Toq — oo corresponds in the numerical calculations to Teq = 10"' fm. 
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Figure 5.1: Time dependence of the asymmetry parameter x (upper part) and the en- 
tropy density normahzed to the values obtained from the Bjorken model a = aoTo/r 
(lower part) for the large initial anisotropy parameter, xq = 100. Thick lines corre- 
spond to the entropy source fl4.22p and thin lines correspond to the source fl5.20p . 
The results are shown for different values of the time-scale parameter T^q- (solid red 
lines), 0.25 fm (dashed green lines), 0.5 fm (dotted blue lines), 1 fm (dashed-dotted 
brown lines), oo (solid black lines). The initial time is fixed to tq = 0.25 fm. 

ii) xo = 0.032(1, where the longitudinal pressure is much larger than the transverse 
pressure, in this case the momentum shape is prolate, and finally, 

iii) the initial system is isotropic, Xq = 1. This case is the closest to the perfect- 
fluid hydrodynamics. In all the cases the system's evolution is studied in the 
time interval 0.25 < r < 10 fm. 

In Fig. I5.1l we show our results for the case i). In the upper part of Fig. I5.1l we show 
the functions x(r) for the entropy source fl4.22p and (15.201) . represented by thick and 

^This value has been chosen because r(lOO) ~ r(0.032) and the case with xq ~ 0.032 may be 
treated as a counterpart of the case xa = 100, with the exchanged roles of the transverse and 
longitudinal pressures. 



42 CHAPTER 5. PURELY LONGITUDINAL BOOST-INVARIANT MOTION 




T [fm] 

Figure 5.2: Same as Fig. 15.11 but the calculations done for the initial anisotropy 

Xq = 1. 

thin lines, respectively. The results are shown for different values of the time-scale 
parameter: Teq = (sohd red lines), Tgq = 0.25 fm (dashed green lines), Xeq = 0.5 fm 
(dotted blue lines), Teq = 1.0 fm (dashed-dotted brown lines), Xgq = oo (solid black 
lines). In the cases Teq = 0.25,0.5, 1.0 fm we observe fast changes of x(r) caused 
by the fact that H{x) and HMsi^) very large for x 3> 1 {H{x) ~ 6x^/^). For 
T ^ Teq the fuuctiou x{t) saturatcs and approaches unity. In addition, we show the 
result for Tgq = 0, which corresponds to sudden isotropization of the system. The 
initial highly-anisotropic system becomes isotropic almost instantaneously. Another 
interesting situation takes place for Teq = oo, where the quadratic behavior (15.291) 
typical for free streaming is observed. In the latter case, the anisotropy grows with 
time and the system never reaches equilibrium. 

In the lower part of Fig. 15.11 we show the corresponding plots of the entropy 
density normalized by the Bjorken solution (15. 7p . During the equilibration process, 
large amount of entropy is produced. For Teq = 0.25, 0.5, 1 fm, the entropy increases 
by about 65%, 80% and 100%, respectively. For r ^ Tgq the normalized entropy 
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saturates, indicating that the flow attains the form of the Bjorken flow. In the case 
where Tgq = oo, we also obtain the Bjorken-like dependence fl5.7p where a ~ A^^/r 
and the transverse temperature becomes independent of time. This is an expected 
result, since without the longitudinal pressure no longitudinal work is done. If, in 
addition, there is no transverse expansion, the system cannot cool down and its 
temperature remains constant. In the = the sudden jump of entropy by 

about 50% is observed. 

In Fig. 15.21 we show our results for the initially isotropic case iii) . The notation 
is the same as in Fig. 15. 1[ In the upper part of Fig. 15. 2[ for Tcq = 0.25, 0.5, 1 fm we 
observe the non-monotonic behaviour of the functions x{t). This is caused by the 
fact that H{x = 0) = and free streaming dominates the early evolution. In this 
case, the anisotropy grows with time {dx/dr > 0). This is an expected feature, since 
even if the starting distribution is isotropic, the system becomes anisotropic if there 
are no interactions to sustain the equilibrium distribution. With increasing r, the 
second therm on the right-hand-side of Eq. (15.51) starts to dominate over the first 
therm and the anisotropy decreases {dx/dr < 0). For r ^ the system becomes 
isotropic and the term 2x/t becomes negligible. In the lower part of Fig. 15.21 we 
observe the corresponding entropy production for Teq = 0.25, 0.5, 1 fm. Again, the 
longer is the anisotropic stage, the more entropy is produced. 

Finally, for completeness, in Fig. l5.3l we show our results for the case ii). For Tgq = 
0.25, 0.5, 1.0 fm we again observe the non-monotonic behaviour of x{t). Initially, the 
dissipation tends to isotropize the system. However, as the equilibrium is reached 
the dissipative processes stop, H{0) = and the free-streaming stage starts to 
dominate the dynamics. The subsequent evolution proceeds very much similarly 
to the case iii). In the lower part of Fig. 15.31 for Tgq = 0.25,0.5, 1 fm, we observe 
the corresponding entropy production. We can see that when the system is near 
equihbrium the entropy production stops. Moreover, contrary to the case i), the 
higher is the value of Tgq, the less entropy is produced. 

Comparing our results presented in Figs. 15. 1[ 15. 2[ and 15.31 we observe that the 
time dependence of x is very much similar in the cases fl4.22p and (15.201) . It suggests 
that the particular form of the entropy source has small effect on the time evolution. 
Far from equilibrium, the cases (I4.22p and (15.201) differ significantly but the large 
entropy production leads to fast decrease of x. As soon as x becomes small these 
two schemes are equivalent, since for small x Eqs. (I4.22p and (I5.20p have the same 
X dependence. 

We also note that Xgq = 0.25,0.5,1.0 fm correspond to t^^ ^ 0.03,0.067,0.13 
fm. Hence, our choice of the parameters corresponds to very short relaxation times 
used in the kinetic theory (in the way done by Martinez and Strickland). This may 
sound reasonable if we consider the plasma as a strongly interacting system. 
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Figure 5.3: Same as Fig. 15.11 and 15.21 but the calculations done for the initial 
anisotropy Xq = 0.032. 



Chapter 6 
Initial conditions 



In order to make realistic comparisons with the data, we should include the trans- 
verse expansion of matter and, possibly, relax the assumption of boost-invariance. 
In this Chapter we introduce the initial conditions used in our model for such situ- 
ations. 

6.1 Glauber initial conditions 
for boost-invariant systems 

In order to solve the ADHYDRO equations discussed in Chapter H] in the boost-invariant 
2-1-1 dimensional [(2-|-l)D] case, it is necessary to specify initial conditions at a 
given transverse space point for the initial proper time r = Tq. The proper 
time r = corresponds to the moment when the centers of the colliding nuclei pass 
through the z = plane. In the numerical calculations we choose Tq = 0.25 fm. 
Before this time, the parton interactions are expected to produce at least a partially 
thermalized system whose later behavior may be described by ADHYDRO. For boost- 
invariant simulations the initial conditions are defined by four functions: cr(ro,Xx), 
a;(ro,x_L), Mx(t-o,x_l), and Mj^(ro,x_L). 

We assume that the initial energy density in the transverse plane, £:o(''"0) ^x)) is 
proportional to the normalized density of sources p(6, x_l), where h is the impact 
parameter corresponding to a given centrality class, 

£o(x±) = £:(ro,x_L) = eip(6,x_L). (6.1) 

The parameter si is the initial energy density at the center of the system in most 
central collisions. Its value may be estimated from the standard hydrodynamic 
calculations. We use the value e\ = 86.76 GeV/fm^. If the considered matter were 
in equilibrium, its temperature would be equal to T = 485 MeV. 

The normalized density of sources, p(6, x_|_), is constructed as a linear combina- 
tion of the wounded-nucleon density pw{b,yi±) and the density of binary collisions 
Pb(&,x±), 

1 — K 

p{b, x_l) = -^—Pw {b, x_l) + kpb (6, xj_) , (6.2) 
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Figure 6.1: Geometry of the collision of two symmetric nuclei. The collision cen- 
trality is specified by the impact vector b. 

The choice of the parameter k {0 < k < 1) controls the mixing between p-yv and 
Pb- In the case k = we obtain the standard wounded-nucleon model [79]. The 
parameter k should be chosen in such a way as to reproduce the experimentally 
observed dependence of the charged particle multiplicity on the collision centrality. 
Typically one uses k = 0.14 for RHIC [HD] and k = 0.15 for the LHC [SI]. 

The wounded-nucleon and binary-collision densities may be obtained from the 
optical limit of the Glauber model. In this approach, for a collision with the cen- 
trality specified by the impact vector b (see Fig. 16.11) . we use the formula for the 
transverse density of wounded nucleons being a sum of the contributions from the 
nuclei with the mass numbers A (projectile) and B (target). Thus, assuming only 
symmetric collisions where A = B and large nuclei {A, B ^ 1) we may write 

Pw ip, x_l) = Pw ^-l) + Pw ^-l) > 

p+(fe,xx) = r(xx-b/2)[l-exp(-ai„r(xi + b/2))], (6.4) 
Pw(&,xx) = r(x^ + b/2)[l-exp(-ai,r(x^-b/2))], 

and the density of binary collisions is 

PB (&, xx) = (J,n T (xx + b/2) r (xx - b/2) . (6.5) 

In Eqs. (16. 4p and (16. 5p Uin is the total inelastic nucleon-nucleon cross section which 
is a function of the collision energy. For the RHIC Au-fAu {A = B = 197) collisions 
at = 200 GeV one chooses a-.^ = 42 mb and for the LHC Pb+Pb {A = 

B = 208) collisions at a/snn = 2760 GeV one can take (Tin = 62 mb [81j . 



6. 1 . GLA UBER INITIAL CONDITIONS FOR BOOST-INVARIANT SYSTEMSA7 




Figure 6.2: Spatial eccentricity ^^(t) = {{y^ — x^))/{{y'^ + x^)) of the wounded- 
nucleon density pw{b) and the binary-coUisions density Pb(&) as functions of the 
impact parameter 6 = |b| calculated for Au+Au collisions at RHIC and for Pb+Pb 
collisions at the LHC 

The role of (Tin is more important in the definition fl6.4p . since it influences the 
shape of profiles. In (16.51) it enters just as a multiplication factor. Thus, the ec- 
centricity of pw{b,^±) for RHIC is larger than that for the LHC, whereas for not 
too peripheral coUisions the spatial eccentricity of pb(^, x^) is in general the same, 
see Fig. 16.21 The density of binary collisions yields a much steeper profile (hard 
component) than the density of the wounded nucleons (soft component). 

The function T(xx) is the nuclear thickness function 

T(x_l) = J dzpws (x±, z) . (6.6) 

The function pws (x±, z) in (16. 6p is the nuclear density function given by the Woods- 
Saxon spherical profile 

Pws (x±, z) = r.^ . ^ T' (6-7) 



1 + exp 



(^^^1 + z^-Ra) /a 



where Ra is the radius of the nucleus with the mass number A calculated from the 
formula 

Ra = 1.12A^/3 - 0.86^-^/^ (6.8) 

a is the surface diffuseness (a = 0.54 fm for both, gold and lead nuclei) and po is 
the nuclear saturation density chosen in such a way that pws (x_l,2;) is normalized 
to the number of nucleons in the nucleus. 

For the initial anisotropy parameter x{tq,:k±), which we take as independent of 
x_L and denote simply as Xq, we consider again three different options: Xq = 100, 
xq = 1, and Xq = 0.032. The case xq = 1 is, of course, the closest to that described 
by the standard perfect-fluid hydrodynamics. The case Xq > 1 corresponds to the 
initial situation where the transverse pressure is much larger than the longitudinal 
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pressure. In this case the momentum shape is oblate (the momentum distribution 
is stretched in the transverse direction and squeezed in the longitudinal direction 
with respect to the beam axis). This type of the initial conditions is considered in 
the Color Glass Condensate (CGC) approach where the distribution functions in 
the longitudinal direction are described by the Dirac delta function, at z = 

[121 SSj, see also Ref. [H]. The values Xq < 1 correspond to the prolate momentum 
shape. This type of the initial conditions has been analyzed, for example, in Refs. 
[581 [83]. 

By fixing both the initial energy density and the initial anisotropy parameter, 
we determine the initial entropy density profile from Eq. 03.341] . namely 



o-o(x±) = o-(ro,x_L) = 



Sipjb, x_l) 
r{xo) 



(6.9) 



Here ^qgpi^) is the inverse function to eqgp{cr). The 6-dependence displayed on the 
right-hand-side of Eq. (16.91) induces centrality dependence of the initial entropy 
density profiles. 

We note that the entropy described by Eq. (16.91) corresponds to the entropy pro- 
duced in the proper time interval < r < tq, just before the ADHYDRO description is 
turned on. For r > Tq, the ADHYDRO model accounts for the additional entropy pro- 
duction described by the entropy source term E. Since the overall entropy produced 
in a given class of collisions is connected with the final multiplicity, we may consider 
it as a fixed quantity. Therefore, if the initial anisotropy xq is much different from 
unity, the initial entropy density (Tq is expected to be small. In this case, the large 
amount of entropy is produced in the phase described by ADHYDRO. On the other 
hand, if the initial anisotropy is close to unity, the subsequent entropy production 
should be negligible. 

Our form of the initial conditions is completed by the assumption that there is 
no transverse flow present at r = tq, i.e., we set Ux(to,x_l) = and Uy(ro,x_L) = 0. 



6.2 Tilted initial conditions 

for non-boost-invariant systems 

The experimental data collected at RHIC show that the boost-invariance is real- 
ized only in a relative narrow central range of rapidities. Thus, the boost-invariant 
(2-|-l)D approach discussed in Section [6711 is not sufficient to perform realistic predic- 
tions for observables outside the midrapidity region (|y| > 1). In order to investigate 
in more detail the influence of partial thermalization on the behavior of matter pro- 
duced in the whole range of rapidities in relativistic heavy-ion collisions, we have 
to perform full-geometry non-boost-invariant 3+1 dimensional [(3+l)D] simulations 
based on the differential equations introduced in Chapter HI 

For a general non-boost-invariant case, the initial conditions for ADHYDRO are 
deflned by flve functions: (T(ro, r], x_l), x(ro, t], x_l), ^^.(ro, r/, x_l), and My(ro, r], x_l). 
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and '(9(ro, r], x_|_). The initial entropy density profile is defined in the similar way as 
in the (2+l)D calculations, 

Mv, x±) = cr(ro, 7], x_l) = ^qgp [^i pib, 7], x_l)] , (6.10) 

however, in the (3+1 )D calculations we choose the parameter Si separately for differ- 
ent physical situations (such as different initial anisotropy profiles and/or different 
values of the time-scale parameter r^q) in such a way as to reproduce well the final 
multiplicity. 

The normalized density of sources p{b,ri,:K±) from fl6.10p . is constructed now as 
a sum of the two terms: 

P(^^,x±) = (1-k) [p+ (6,Xx)/+(r/)+p^(6,Xx)/-(r/)] 

+ «:pB(&,xx)/(r7). (6.11) 

Here, the function / (77) is the initial longitudinal profile 

{rj — Arjy 



f [r]) = exp 



2(t2 



-o{\v\-^v) 



(6.12) 



The first term in fl6.1ip introduces contributions to the density of sources from the 
forward- (-I-) or backward-moving (— ) participant nucleons. This formula assumes 
a preferred emission from the participating nucleons along the direction of their 
motion. The second term introduces a symmetric contribution from the binary 
collisions. The half-width of the central plateau. At], and the half-width of Gaus- 
sian tails, fitted to reproduce the RHIC multiplicity of charged particles in 
pseudorapidity. 

The initial density of sources produced by a single forward- or backward-moving 
participant nucleon with rapidity yh = ln(y'sNN/'"^N), where rriN is the nucleon mass, 
is proportional to 

r^'Hv) = f iv) fF(B) iv) , (6.13) 

where 

rj < -rjr, 

MV) = { ^ -V.n <V<Vm, fsiv) = fri-v)- (6.14) 

,1 Tlm<1l 

The range of rapidity correlations equals = y;, — ?7s ^ 3.36 and the shift in the 
rapidity 77^ = 2 is a phenomenological parameter. With the assumption fl6.14p . the 
hydrodynamic description outside the region {—rjm, rjm) is not rehable. 

Similarly to (2+l)D calculations, for the initial anisotropy profile, x(ro, r], x^), 
we choose three different values: Xq = 100, xq = 1, and Xq = 0.032. Moreover we 
check the possibility of the spatial dependence of the initial anisotropy parameter, 

= Xoi.pi.Vy x±))- We assume also that there is no transverse flow present initially, 
therefore, we set Ma(To, ''7, x_l) = and Mj,(ro, ?], x_l) = 0, and the initial longitudinal 
rapidity of the fluid follows the simple Bjorken scaling '(9(ro, r/, x^) = rj. 
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Chapter 7 



Freeze— out 



In this Chapter we present first the method used in ADHYDRO to extract the freeze-out 
hypersurface from hydrodynamic evolution. Then, the Cooper-Frye prescription in 
the form used by the Monte-Carlo event generator THERMINATOR [Ml [55] to describe 
particle production on a given freeze-out hypersurface is derived. Finally, calcula- 
tions of physical observables with the help of THERMINATOR are explained in more 
detail. 

7.1 Single- freeze-out approximation 

The rapid expansion of dense and hot matter (described in our approach by ADHYDRO) 
implies that the density of particles decreases and the mean free path of particles 
grows. When the time between interactions of particles becomes larger than the 
expansion rate (described by divergence of the four-velocity field) 



the collision processes cannot maintain the local equilibrium any longer. The strongly 
interacting system of particles undergoes a transition to the system consisting of es- 
sentially free-streaming particles. This complex dynamic process, called freeze- out, is 
quite often treated approximately as if it took place locally on the three-dimensional 
freeze-out hypersurface. 

Typically, the freeze-out condition (17. ip is realized by assuming that hadrons 
completely decouple at the space-time points x where the entropy density drops be- 
low a certain freeze-out entropy density cxf. For systems in local thermal equilibrium 
described by the equation of state o"(T) this condition may be formulated in terms 
of temperature. 



A set of freeze-out points x extracted from the hydrodynamic evolution by the 
constraint (17. 2p forms the freeze-out hypersurface. One can subsequently use it to 
generate primordial particles employing the Cooper-Frye formalism. 



Tcoll > 




T{a{x)) = T, = T{a,). 



(7.2) 
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In general, the inelastic processes in the expanding system cease earlier than 
the elastic ones. Therefore, one can distinguish chemical freeze-out (the stage when 
inelastic processes stop and hadron abundances are fixed) from subsequent kinetic 
freeze-out (the stage when all types of interactions cease and the momentum dis- 
tributions of hadrons are frozen). Of course, the latter may be considered as the 
true freeze-out process. In the case of high energy collisions at RHIC, which are 
the subject of our studies, we can use quite efficient simplification called the single- 
freeze-out approach [861 EZl EBJ. In this approach, one assumes that the kinetic 
freeze-out coincides with the chemical freeze-out, Tkin = Tchem = Tf. In our studies 
we adopt this scenario. 



7.2 Cooper- Frye formalism 

The momentum distribution of primordial hadrons emitted on the freeze-out hyper- 
surface S may be calculated from the Cooper-Frye formula [89] 



dN 



dyd'^PT 



g J dE,{x)p^f{x,p), (7.3) 



where dE'^ is the element of the three-dimensional freeze-out hypersurface submerged 
in the four-dimensional Minkowski space, f{x,p) is the phase-space distribution 
function of the emitted particles and ^ = 2s + 1 is the spin degeneracy factor, where 
s denotes the particle's spin. For systems in local thermal equilibrium described by 
the perfect fiuid hydrodynamics one assumes equilibrium Fermi-Dirac (e = —1) or 
Bose-Einstein (e = +1) statistics, i.e., one uses the distribution functions 



f{x,p) = < exp 



p^,U^{x) - /i' 



(7.4) 



depending on the scalar product of particle four-momentum p^^ and the four- velocity 
If^ of the fiuid element. Clearly, if deviations from equilibrium at freeze-out are 
large, the distributions (17.41) should be replaced by the appropriate non-equilibrium 
distributions. In our calculations, for a typical choice of parameters, only initial 
parts of freeze-out hypersurfaces exhibit large anisotropics. Thus, as long as they 
weakly contribute to total final multiplicities they may be simply neglected in the 
calculation of the observables. 

One can determine thermodynamic parameters on the freeze-out hypersurface 
with the help of thermal models [90l [9T| [921 [93l IM] which analyze relative hadronic 
abundances. At RHIC, this type of study [91] indicates the non-zero values of 
chemical potentials. Thus in the generation of hadrons we introduce the chemical 
potential fi 

/i = BfiB + hfJ'h + Sfis + Cfic, (7.5) 

where B, J3, S, C is the baryon number, third component of isospin, strangeness, 
and charm of the emitted hadron, respectively. Similarly, /i^, /i/3, /i5, and nc denote 
the corresponding chemical potentials. 
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The particle's four-momentum may be parametrized in the standard way 

= {rriT cosh y, px cos 0p, px sin 0p, rriT sinh y) , (7.6) 

where mr = + p^ is the transverse mass, pr is the transverse momentum, y is 
the (longitudinal) rapidity of the emitted particle, and (f)p is the azimuthal angle of 
the transverse momentum. Together with the parametrization (14.51) . one may write 



P^lU^ = ^1 + ul + ul rriT cosh(^9 - y) - Pt{ux cos 0p + Uy sin 0p). (7.7) 

The form of dS^ may be obtained with the help of the formula known from differ- 
ential geometry [95] 

dT.^ = Ef^ais-f -g^ —dadhdc, (7.8) 

where e^ais-) is the antisymmetric Levi-Civita tensor with ^0123 = +1- The quantity 
dTi^{x) defines a four- vector that is perpendicular to the hypersurface at point x. 
Its norm is equal to the volume of the hypersurface element. The variables a, h and 
c introduce a coordinate system in Minkowski space parameterizing the positions of 
points on the freeze-out hypersurface. Their ordering is chosen in such a way that 
dE'^ points in the direction of decreasing temperature. 

The freeze-out hypersurface obtained from the most general (3+l)D hydrody- 
namic evolution (without assuming any symmetries) may be parametrized in the 
following way [11] (see Fig. 17.11) 

, I n\ ■ n ■ ^\ 1 d(C,(l),9) cos9 
t = (tq + a(C, 0, t^) smfc'smC) cosh — , 

X = d{(, (f),9) sinO cos( cos(f), 
y = d{(, (f),9) sinO cos(sm(j), 

z = (ro + ci(C,0,^) sin ^ sin sinh ^^^1^^^^^^. (7.9) 

This parametrization leads to simple formulas for the space-time rapidity rj, longi- 
tudinal proper time r, and the transverse distance r 0, 

d(C,(f),e) cosO 

^ = A ' 

r = tq + d{(, (j),9) sinO sin(, 

r = rf(C,0,^)cosC. (7.10) 

The parameters a, h and c in Eq. (17. 8p are identified with three angles and 9, 
respectively, which are restricted to the ranges: 

< C < 7r/2, 
< < 27r, 

< 6 < TT. (7.11) 



^In the calculations we remove the emission points with r < 1 fm/c. In this region the system 
is anisotropic and consists mainly from gluons. 
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Figure 7.1: Coordinate system used to parametrize the freeze-out hypersurface extracted 
from a (3+l)D dimensional hydrodynamic code, see Eqs. (|7.9p - (|7.11|) . The part (a) shows 
the view in x — y — t space with the fixed value of 9, while the part (b) shows the view in 
r — Arj — T space with the fixed value of (p. 



The quantity d{(, 0, 6) describes the distance between a point on the freeze-out hy- 
persurface and the coordinate system's origin {t = tq,x = 0,y = 0,1] = 0). The 
parameter Tq is the initial proper time for hydrodynamic evolution, at which the 
initial conditions are specified. The parameter A introduces a scale which changes 
dimensionless space-time rapidity r] into a quantity with dimension of fermi. The 
parametrization (17. 9p works for all typical freeze-out hypersurfaces where the dis- 
tance d is a function of (, (j), and 9. 



Substitution of (17. 9p in (17.81) leads to the following form of the Cooper-Frye 
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integration measure 
(fr sin 9 



A 



dd . , . ^ 
— cos Q [Pt sm Q cos I 



dd 



,) — rriT cos C cosh — r])) 



+ cos C sin 6* ( (i sin — — cos 9 j {pt cos C cos (0 — 0p) + sin ^ cosh — rj)) 

+ cos C sin 6* ( d cos 6* + ^ sin 9] —rriT sinh — rf) 

89 J T 



(^d . , , , , 
+ ^ Pt sin(0 - 0p) 



dCd(j)d9. 



(7.12) 



Using formulas (17.71) and (I7.12p in (17.31) we obtain obtain a six-dimensional parti- 
cle distribution which one can use to generate both stable hadrons and unstable 
resonances on the freeze-out hypersurface. 



7.3 Extraction of physical observables 

After performing hydrodynamic calculations with ADHYDRO, the freeze-out hyper- 
surface is extracted. This procedure yields four functions: d{(,(f),9), U:i:{(,(f),9), 
Uy{(,(f),9) and ^{(,(f),9) which are subsequently used by the THERMINATOR Monte- 
Carlo generator as an input. THERMINATOR generates physical events using the 
Cooper- Frye formula (17.31) . In the single event all known hadron species (stable 
and unstable) are generated on the freeze-out hypersurface. Subsequently, the un- 
stable primordial particles decay through strong and electro-weak interactions. The 
final particle distributions consist of stable primordial particles and contributions 
from resonance decays. The complete information about all generated particles and 
their decays is kept in the program. This information allows for the calculation of 
soft hadronic observables typically measured in heavy-ion experiments. 

Probably, the most important information about the evolution of matter created 
in relativistic heavy-ion collisions is contained in the momentum distributions of 
identified particles. Experimental data on transverse-momentum distributions are 
often corrected for contributions from weak decays (the contributions from weak 
decays to the hadronic yields and spectra are experimentally determined and sub- 
tracted). The detailed output of THERMINATOR allows for applying such corrections 
in the model calculations in the easy way, making comparisons to the data more 
realistic. 

In non-central collisions of two heavy nuclei, the experimentally observed mo- 
mentum distribution of particles is azimuthally anisotropic and may be expressed as 
a Fourier series in the momentum azimuthal angle, (pp = SiTctsin[py/px), with respect 
to reaction plane 



d^N 



d^N 



d^Pxdy 2npT dpxdy 



1 + ^ ^VnipT, y) cos{n(j)p) 



(7.13) 



n=l 
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In the hydrodynamic calculations where the initial conditions are defined by a 
smooth density profile taken from the optical Glauber formula fl6.2l) the sine terms 
disappear due to symmetry with respect to the reaction plane. In event-by-event 
calculations where the fiuctuations of initial conditions are present this is not the 
case. The coefficient Vn{pT,y) is the n-th harmonic differential flow which quantita- 
tively characterizes the momentum azimuthal anisotropy. The Vi{pt,y) coefficient 
is called the directed flow, and the V2{pt, y) coefficient is called the elliptic flow. Due 
to symmetry reasons the vi{pt, j) vanishes for y ~ 0. The coefficient is called the 
triangular flow. The higher harmonics are typically negligible. For instance, the f 4 
coefficient is predicted by hydrodynamic calculations to be at the level of ~ 0.1%. 

Having complete information about the particles generated by the THERMINATOR 
code one can perform the femtoscopic analysis. The calculation of the correlation 
function for the identical pions is performed using the two-particle Monte-Carlo 
method described in Ref. [HZl EH] and implemented in THERMINATOR 2 [85] code. In 
the numerical procedure the ideal Bose-Einstein correlation function is evaluated 
using the formula 

EE'^A(q-p. + P,) 5A(k-i[P. + P.]) I^(q*,r*)|' 

= ^ '''VvA ( ^ w M ^ n • ^^-^^^ 

L L (q - Pi + Pj) ()a (k - 2 [Pi + Pj\) 

i j^i 

The calculation of the correlation function (17.141) is performed with the help of the 
box function (5a defined in the following way 

lO, otherwise, 

where one uses the bin resolution A = 5 MeV. In Eq. (I7.14p we use the momentum 
difference of the pions, 

g = (?o, q) = (^pi - Ep^,Pi - P2) , (7.16) 
and the average momentum of the pion pair, 

k = {ko, k) = ^ {Ep, + Ep„p, + P2) . (7.17) 
We also consider the space-time separation between pions 

Ax = (t,r) = (ti-t2,ri-r2). (7.18) 



In the HBT analysis we use the standard Bertsch-Pratt system of coordinates 
1100] ■ where the long direction is parallel to beam axis, the out axis coincides with 
the transverse momentum of the pion pair kx and the side axis is perpendicular to 
out and long axes. We perform the analysis in the longitudinally co-moving system 
(LCMS) in which fciong = 0. For each pair of pions the wave function in Eq. fl7.14p . 

^{ci*,r*) = ^{e^^*^'' +e-'^*^*), (7.19) 
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is calculated in the pair rest frame (PRF). The q* and r* denote the relative momen- 
tum and position calculated in PRF. In the calculation of expression (17. 141) we do 
not include Coulomb effects. Using the Bertsch-Pratt coordinates kr, (?out, Qside, Q'long 
the correlation function (17.141) may be approximated by the formula 

C {k±, gout, gside, giong) = 1 + A CXp [-i?out )?out 

-RLe{k±)qLe - RL,ik±)qL,] , (7-20) 

where the HBT correlation radii -Rout, -Rside and -Riong are the widths of the Gaussian 
source. From the fit we obtain parameters -Rout, -Rside and -Riong as a function of kx- 
In order to make a comparison of our results with the available data from STAR we 
consider four kr bins (0.15-0.25), (0.25-0.35), (0.35-0.45), and (0.45-0.60) GeV. 
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Boost-invariant systems with 
transverse expansion 

The data from RHIC and the LHC shows that the observables measured in the cen- 
tral rapidity region are approximately invariant with respect to the Lorentz boost 
transformations along the beam axis. Thus, in this Chapter, we use ADHYDRO as- 
suming the boost-invariance along the beam direction. In this way, we restrict our 
analysis to the midrapidity region (|y| < 1) of ultra-relativistic heavy-ion collisions. 
We perform calculations for central, c = — 5% (6 = 2.26 fm), and mid-central col- 
hsions, c = 20 - 30% {b = 7.16 fm) and c = 20 - 40% {b = 7.84 fm), and compare 
the model results with the available data from Au-|-Au collisions at the energy of 
y^SNN = 200 GcV at RHIC. We check whether the existence of high anisotropics 
of pressure at the early stages of evolution may be probed with the midrapidity 
observables. 

8.1 Central collisions 

First we consider central Au-|-Au collisions at the energy of ^snn = 200 GeV 
at RHIC. We choose the ccntrality class c — — 5%, which corresponds to the 
impact parameter b = 2.26 fm. We check three scenarios of the initial momentum 
anisotropy: 

i) xq = 100 (dashed blue lines), which means that the system is initially oblate in 
the momentum space, the fact predicted by the microscopic models of heavy- 
ion collisions, 

ii) Xq = 1 (solid black lines), which corresponds to initially isotropic systems in 
the momentum space, this case gives the closest description to perfect fluid 
hydrodynamics, and 

iii) Xq = 0.032 (dotted green lines), which means that the system is initially prolate 
in momentum space, which is less-likely to be realized in nature than the case 
i), however, we consider this case for completeness. 
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Figure 8.1: Left panel: Time dependence of the central entropy density for the centrality 
class c = — 5%. Results obtained for three values of the initial anisotropy parameter: 
xq = 100 (dashed blue line), xq = 1 (solid black line), and xq = 0.032 (dotted green line) as 
discussed in the text. The results are compared to the reference perfect-fluid calculations 
(dashed-dotted red line). Right panel: Time dependence of the entropy production in the 
center of the system. The prefect-fluid hydrodynamics is entropy conserving, hence, it is 
not included in the plot. 



The results of the ADHYDRO model are compared to the reference perfect-fluid cal- 
culations performed with the LHYQUID code [1011 182] . In our calculations we use 
the initial conditions described in Chapter 16. 1[ The value of the initial energy 
density Si = 86.76 GeV/fm^ in Eq. (16.11) is the same for all analyzed cases. For 
the system in local thermal equilibrium, it corresponds to the initial central tem- 
perature Ti = 485 MeV. Following Eq. (16. 9p . the initial central entropy density is 
(Ti = cr(ro, = 0) = 227 fm~'^ for the case xq = 1, and cXi = 151 fm~^ for the cases 
Xo = 100 and xq = 0.032. The value of the time-scale parameter Tcq = 0.25 fm used 
in our calculations guarantees that the system reaches local thermal equilibrium in 
the whole fireball volume approximately after 1 fm (this value for thermalization 
time is easily acceptable by most of the microscopic model calculations). 

In the left panel of Fig. l8.1l we present the time dependence of the central entropy 
density a{T, x± = 0) for three considered cases of the initial anisotropy parameter: 
Xo = 100 (dashed blue line), Xq = 1 (solid black line), and Xq = 0.032 (dotted green 
line). We compare our results to the perfect-fluid case (dashed-dotted red line). 
We observe that in the cases xq = 100 and xq = 0.032 the entropy density grows 
initially. This behavior is different from the entropy-conserving perfect-fluid case 
where (T(r, x_l = 0) decreases like 1/r (dashed-dotted red line). For r ^ 0.3 fm 
(j(r, xx = 0) reaches maximum and for r > 0.5 fm it starts to behave similarly as 
in the perfect-fiuid case. The case xq = 1 gives results similar to the perfect-fiuid 
case. The entropy densities at the end of evolution are similar, however, small 
differences are noticeable. It means that the rescaling of the initial entropy density 
(Ti according to Eq. (16. 9p assures only approximately the same final multiplicity and 
entropy density in all considered cases. 
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T [fm] T [fm] 

Figure 8.2: Left panel: Time dependence of the anisotropy parameter at the center of 
the system for the centrahty class c = — 5% for three values of xq. Notation is the same 
as in Fig. 18. 1[ Right panel: In-plane ((/> = 0) time dependence of the radial velocity v± for 
different values of radial distance. 

In the right panel of Fig. 18. II we show the time dependence of the entropy source 
E(r, x_|_ = 0) defined by Eq. fl4.22p . The notation is the same as in the left panel. 
The function E(r, = 0) contains information complementary to that encoded 
in (j(r, = 0). Due to strong initial anisotropy, in the case xq = 100 the initial 
entropy production is very large. As the pressure anisotropy decreases (see the 
corresponding time dependence of the anisotropy parameter x{t, x^ = 0) in the left 
panel of Fig. l8.2p the production of entropy becomes smaller. In the case Xq = 0.032, 
the entropy production is also quite large initially. However, it quickly drops to 
zero when x passes unity (see also Chapter 15.61) . Subsequently, it increases again 
and becomes similar to the values found for the case xq = 100 0. When xq = 1, 
the anisotropy remains always close to unity. Hence, the entropy production is 
very small in this case. We note that in all the cases the anisotropy parameter 
tends asymptotically to unity, which means that for large times the entropy source 
vanishes and the system approaches the perfect-fiuid regime. 

The initial large anisotropy in the case Xq = 100 is equivalent to the large P±/P\\ 
ratio which implies strong hydrodynamic flow in the transverse direction at the early 
stages of the evolution (small longitudinal pressure means that no work is done in 
the longitudinal direction and, consequently, thermal energy is used to build the 
transverse flow only). This effect is, however, compensated by the smaller value of 
the initial entropy density obtained from Eq. (16. 9p . Thus, our results obtained with 
ADHYDRO are very much similar to those obtained in the standard hydrodynamics. 
This property is seen in the right panel of Fig. 18. 2^ where we plot in-plane (</> = 0) 
radial velocity v± as a function of proper time r for different values of the distance 
r, for the cases i), ii), and iii) H. 



^Different entropy productions in the cases xq — 100 and xq — 0.032 are responsible for slightly 
different final entropy densities. 

^Since for c — — 5% the evolution is almost azimuthally symmetric in the transverse plane, 
the out-of-planc {(f) ~ 7r/2) behavior is very similar to the one shown in the right panel of Fig. 18.21 
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Figure 8.3: In-plane (0 = 0) and out-of-plane (0 = it/2) freeze-out hypersurfaces for the 
cases xq = 100 (dashed blue lines) and xq = 0.032 (dotted green lines) shown together 
with the values of the velocity w±(C) for the centrality c = — 5%. ADHYDRO results are 
compared to the perfect-fluid case (dashed-dotted red lines). The freeze-out temperature 
equals Tf = 150 MeV. 



In Fig. 18.31 we present in- and out-of-plane freeze-out hypersurfaces determined 
from the numerical solutions of the ADHYDRO equations. The freeze-out condition is 
specified only for the locally equilibrated matter, thus, we exclude the space-time 
regions where r < 1 fm. In those regions the matter is anisotropic and consists 
mainly of gluons — it is not likely that it contributes significantly to the hadron 
production. 

We use the freeze-out temperature value Tf = 150 MeV which has turned out 
very much satisfactory to reproduce the experimental data from RHIC using perfect- 
fluid hydrodynamics. After extracting the freeze-out hypersurface from the hydro- 
dynamic code, we use THERMINATOR 2 [85] in order to generate physical events. In 
Fig. l8.3l we can see that for xq = 100 the system cools down a little bit slower than in 
the standard perfect-fluid hydrodynamics, so the freeze-out hypersurface is larger. 
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Pt [GeV] Pt [GeV] 

Figure 8.4: Left panels: The calculated transverse-momentum spectra of vr"*" (upper part), 
(middle part), and protons (lower part) for Au+Au collisions at -^snn = 200 GeV 
and the centrality class c = — 5% compared to the PHENIX data |102j . Protons are 
feed-down corrected for A°. Right panels: The calculated transverse- momentum spectra 
of hyperons A'' (upper part), H~ (middle part), and + Cl~ (lower part) for Au-|-Au 
collisions at y^SNN = 200 GeV and the centrality class c = — 5% compared to the STAR 
data [103j . A'' corrected for feed-down from H and $7. 

At the same time we find that the values of the hydrodynamic flow are slightly 
larger. This, as will be shown later, leads to a slightly larger multiplicity and harder 
spectra. In the case Xq = 0.032, we deal with the opposite situation. Despite these 
small differences, in all the cases the freeze-out hypersurfaces are similar. This sug- 
gests that the initial anisotropic stage of evolution does not modify significantly the 
freeze-out conditions, provided the initial entropy density is appropriately rescaled 
in each case. 

In the left panels of Fig. 18.41 we show the transverse- momentum spectra of pions 
(upper part), kaons (middle part), and protons (lower part), while in the right 
panels of Fig. 18.41 we present the transverse-momentum spectra of A° (upper part), 
S~ (middle part), and fl~^ + fl~ (lower part). One can observe that for all values 
of the initial anisotropy the spectra are very similar and lie close to the (2-|-l)D 
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Fi gure 8.5: HBT radii: i?out5 -^sidei -^long ^ ^ function of total transverse momentum 
kT of the pion pair, calculated for positive pions with THERMINATOR 2 [85j for freeze-out 
hypersurfaces coming from ADHYDRO and LHYQUID. The results are presented together with 
the STAR data (dots). 



perfect-fluid hydrodynamics results. The small differences can be noticed in the 
high pt region where the results for xo = 100 are the highest and for xq = 0.032 
are the lowest. This is mainly due to larger values of transverse velocity of the fluid 
at the freeze-out (see discussion of Fig. I8.3p and slightly larger multiplicity (see left 
part of Fig. 18. ip in the case (i). 

We close the study of the central collisions with Fig. 18.51 which presents the 
results for /cT-dependent HBT radii. We observe the overall agreement with the 
experimental data for the perfect-fluid calculation performed with LHYQUID. We 
also see that the correlation radii show almost no sensitivity to the initial anisotropic 
stage. All calculations yield similar results. This is understandable, since we obtain 
similar shapes of the freeze-out hypersurfaces for all considered initial values of Xq. 

The analysis of the most central collisions in the framework of ADHYDRO and 
LHYQUID suggests that the considered observables are almost insensitive to the initial 
anisotropic stage as long as the initial entropy density is properly rescaled. 

8.2 Non-central collisions 

In this Section we turn to the analysis of non-central collisions. We take into consid- 
eration two centrality bins: c = 20 - 30% (b = 7.16 fm) and c = 20 - 40% (b = 7.84 
fm). The same value of the initial central energy density, ei = 86.76 GeV/fm^, is 
used now, as before for central collisions. The initial entropy density at the center of 
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Figure 8.6: Same as Fig. [O but for the centrality class c = 20 - 30% (6 = 7.16 fm). 

the system is rescaled according to Eq. (16. 9p . In this way we find: a-i = 183.1 fm^^ 
for the case Xq = 1, and ai = 120.3 fm~^ for the cases Xq = 100 and Xq = 0.032 if 
c = 20 — 30% [similarly: 0"; = 171.3 fm~^ for the case Xq = 1, and o"i = 113.8 fm^^ 
for the cases Xq = 100 and Xq = 0.032 if c = 20 - 40%]. 

In the left part of Fig. 18.61 we present the in-plane (0 = 0) and out-of plane 
(0 = 7r/2) freeze-out hypersurfaces for the centrality class c = 20 — 30%. We come 
to similar conclusions as in the case of central collisions. Moreover, we observe that 
the pressure acts in the same way in the in-plane and out-of-plane directions. This 
feature is expected since we introduced the anisotropy only between the longitudinal 
and transverse pressures. We also observe that the transverse spatial eccentricity of 
the fireball, e^, decreases faster with time in the Xq = 100 case than in the Xq = 0.032 
case, see Fig. 18 .71 This implies that the transverse momentum anisotropy Sp gets 
built faster. However, these effects are small if the equilibration time is short. 

In the right panel of Fig. 18.71 we show the differential elliptic flow coefficient 
V2{Pt) of pions as a function of the transverse momentum, for the centrality class 
c = 20 — 30%. The model results describe quite well the data for Pt < 1 GeV. 
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Figure 8.7: Left panel: The spatial eccentricity ^^(t) = {{y"^ — x^))/{{y^ where the 

average is done with the energy density as a weight function, and the transverse-momentum 
anisotropy Spij) = (/ dxdy{T^^ — Tyy))/{J dxdy{T^^ + T^^)) shown as functions of time. 
Calculation are done for centrality class c = 20 — 30%. Right panel: The elliptic flow of 
pions in Au+Au collisions at -^/snn = 200 GeV and the centrality class c = 20 — 30%. The 
model results are compared to the RHIC experimental data |105] . 

For different initial anisotropy parameters we obtain almost the same results. The 
slightly larger values of V2 are observed in the case xq = 100, which is consistent 
with the conclusions coming from the analysis of the parameter Sp. However, for 
small Teq that is used by us, this effect is negligible. 

In Fig. 18.81 we show our model results for the transverse-momentum spectra of 
pions (upper left part), kaons (middle left part), protons (lower left part), A° (upper 
right part), S~ (middle right part), and + (lower right part) obtained for the 
centrality classes c = 20 — 30% and c = 20 — 40% and compared with the RHIC 
data from PHENIX [1Q2J and STAR |103] . One can notice the very good agreement 
of the model results with the experimental data, except for A*^ whose normalization 
is too large. The model describes the experimental data up to px = 2 GeV, thus 
it works well in the region of applicability of hydrodynamics. The discrepancies 
for large pt are more significant than in the central collisions. Similarly to central 
collisions we observe improvement of the high pt tails for the case of xq = 100. 

Figure 18.91 presents the dependence of the HBT radii for the centrality class 
c = 20 — 30%. Similarly to central collisions the agreement with the data is quite 
satisfactory. There is almost no effects of the early anisotropic stages. 

Summarizing, the systematic study of the observables in the (2+l)D boost- 
invariant version of ADHYDRO model brings us to the conclusion that the observables 
are almost insensitive to the initial highly-anisotropic stage, provided the initial 
entropy density is properly rescaled. Such results are consistent with the universality 
of the flow predicted in Ref. |106] . 
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Figure 8.8: Left panels: The calculated transverse-momentum spectra of vr"*" (upper part), 
(middle part), and protons (lower part) for Au+Au collisions at -^/snn = 200 GeV and 
the centrality class c = 20 — 30% compared to the PHENIX data (dots) [102J. Protons are 
feed-down corrected for A''. Right panels: The calculated transverse- momentum spectra 
of hyperons A*^ (upper part), H~ (middle part), and + il^ (lower part) for Au+Au 
collisions at y^SNN = 200 GeV and the centrality class c = 20 — 40% compared to the 
STAR data (dots) [103]. A'^ corrected for feed-down from H and fi. 
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Figure 8.9: HBT radii: iiout, Rside, -Rlong, calculated with ADHYDRO and LHYQUID. 
The model results for Au+Au collisions at -^/snn = 200 GeV and the centrality class 
c = 20 — 30% are compared to the STAR experimental data (dots) taken from |104j . 



Chapter 9 



Non-boost-invariant description 
of RHIC data 



The RHIC results obtained outside the midrapidity region [Wl [IMl IH IIDHI fTTUl 
1111] . in particular, the measured rapidity distributions of charged particles, suggest 
that any realistic description of the data should relax the assumption of Bjorken 
boost-invariance. Following this observation, in this Chapter we use the ADHYDRO 
model in 3+1 dimensions. We do not impose any symmetries to analyze the effects 
of initial anisotropic stages on the subsequent behaviour of matter. The analysis 
is done for Au+Au collisions studied at RHIC at the energy ^snn = 200 GeV. 
In order to make a comparison with the experiment, we use the ADHYDRO model 
coupled to the Monte-Carlo model of statistical hadronization THERMINATOR [Ml 185] 
as described in Chapter [71 

9.1 Constant initial momentum anisotropy 

First we analyze the simplest case where the initial anisotropy Xq is constant in 
space. Once again, we perform calculations for three different cases: i) xq = 100 
(dashed blue lines), ii) xq = 1 (solid black lines), and iii) xq = 0.032 (dotted green 
lines). We keep this notation all the way in this Chapter. We also use three values 
of the time-scale parameter: Tgq = 0.25 fm, Tgq = 0.5 fm, and Tgq = 1.0 fm to check 
how the results differ with respect to the length of the anisotropic stage. We show 
the results for the fixed freeze-out entropy density Uf = 1.79/fm^ which corresponds 
to the freeze-out temperature Tf = 150 MeV. This value of Tf correctly describes 
the slope of pr spectra of different hadron species, thus satisfactory separating the 
flow component from the thermal component in transverse- momentum spectra. The 
half-width of the central plateau in the initial profile, Ai] = 1, and the half-width 
of gaussian tails, = 1.3, are fixed, unless specified otherwise. They are fitted to 
approximately reproduce the shape of charged particle distributions in pseudora- 
pidity in the case xq = 1. The initial energy density in the center of the fireball, 
Ei, is chosen in such a way as to reproduce the total experimental charged particle 
multiplicity A^^h'' ~ 5060 ±250 [107j for the centrality class c = — 6%. For the case 
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Xo 1 100 0.032 

Teq [fm] 0.25 0.25 0.5 1.0 0.25 0.5 1.0 

£i [GeV/fm3] 107.5 58.3 53.4 48.8 72.5 75.9 80.1 



Table 9.1: Values of the initial central energy density, e^, used in all the figures of this 
Chapter except for Fig. 19.21 



Xo = 1 we use the energy density Ei = 107.5 GeV/fm^ which yields N^^'^°'^ = 5020. 

Since entropy grows during the hydrodynamic evolution due to non-zero entropy 
production source S(req, a, x), see Eq. (14.221) . the initial central energy density must 
be appropriately renormalized for each pair of parameters x and Teq to keep N^^"^"^ 
unchanged. In the case Xq = 100 the total entropy produced during the hydrody- 
namic evolution grows with increasing Tgq, while in the case Xq = 0.032 the total 
produced entropy decreases with increasing Teq (similarly to the one-dimensional 
behavior shown in Figs. 15.11 and 15.31 in Chap. 15.61) . 

The values of ei are summarized in Table 19.11 The centrality dependence is 
reproduced by applying the scaling (16.31) . The initial conditions are given by the 
tilted source described in Chapter 16. 2[ 



9.1.1 Particle spectra 

We start our discussion of the results obtained with the ADHYDRO model with the 
analysis of the momentum distributions of charged and identified particles. These 
distributions include the feeding from unstable resonances whose decays are imple- 
mented in THERMINATOR [Ml 185] . 



Pseudorapidity distribution 

At first, we consider the pseudorapidity distributions of charged particles. The 
pseudorapidity is defined by the formula 

rjps = -In , (9.1) 

2 p-pii 

where p is the momentum of a particle. In Fig. 19.11 we show dNch/drjps for three 
different values of the initial momentum anisotropy: Xq = 100 (dashed blue lines), 
Xo = 1 (solid black lines), and xq = 0.032 (dotted green hues), and three values of 
the time-scale parameter: Teq = 0.25 fm, Teq = 0.5 fm, and Teq = 1.0 fm. We observe 
that as long as the initial anisotropic stage is short its influence on pseudorapidity 
distributions of particles is small in both xq = 100 and Xq = 0.032 cases, as com- 
pared to the case Xq = 1. As the length of the off-equilibrium stage increases, the 
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Figure 9.1: The pseudorapidity distribution of charged particles for three values of the 
time-scale parameter: Teq = 0.25 fm (a), Teq = 0.5 fm (b), Teq = 1.0 fm (c). The results 
are obtained for three values of the initial momentum anisotropy: xq = 100 (dashed blue 
lines), xo = 1 (solid black lines) and xq = 0.032 (dotted green lines) and two centrality 
classes: c = — 6% (6 = 2.48 fm) and c = 15 — 25% (6 = 6.4 fm). The results are compared 
to experimental data from PHOBOS |107j (red dots) and BRAHMS [108j (green squares). 
Vertical bars for PHOBOS data denote the systematic errors. 

difference of pressures lasts longer, and the difference between the longitudinal and 
transverse expansion rates becomes larger. For instance, in the case xq = 100 the 
pseudorapidity distribution becomes a bit steeper in the region of forward and back- 
ward rapidities due to reduced longitudinal pressure (see our discussion in Chapter 
15. 3p . but the width of the central plateau is nearly unchanged. 

In Fig. 19.21 we present the case Xq = 1 with Tpq = 0.25 fm and the standard pa- 
rameters of f{r]) together with the case Xq = 100 with Teq = 1.0 fm. In the latter case 
the initial energy density was reduced to = 41.8 GeV/fm'^ (from 48.8 GeV/fm^, 
see Table W/Lt) and the width of the plateau of the initial profile f{ri) was refitted to 
Arj = 1.5. In this way we obtain the same total multiplicity. We can clearly see 
the improvement in the description of data as compared to the case Xq = 100 in 
Fig. 19.11 (c). Thus, the shape of pseudorapidity distributions might be better de- 
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Figure 9.2: The pseudorapidity distribution of charged particles for xq = 1 (soUd black 
line) and xq = 100 (dashed blue line). In the latter case the initial conditions were modified 
to Si = 41.8 GeV/fm'^ (compare Table [9?T]) and Ar/ = 1.5 to obtain the same final total 
multiplicity. The results are compared to the experimental PHOBOS |107j (red dots) and 
BRAHMS fLU8\ (green squares) data. 

scribed if the initial fireball is transversally thermalized and the profile parameters 
are properly readjusted. However, it should be stressed that the widths and forms 
of the initial longitudinal profiles fl6.12p are not well known. Therefore, no definite 
conclusions about the length of the pre-equilibrium phase may be drawn from the 
correct description of the shape of the rapidity distributions alone, as it has been 
already noticed in [112] . 

Transverse-momentum spectra 

In Fig. 19.31 we present the transverse- momentum spectra of positive pions (left part) 
and kaons (right part) for the centrality class c = — 5% (6 = 2.23 fm) and for 
different rapidity windows. The model calculations are done with Teq = 1.0 fm. The 
model results are compared to the BRAHMS data [Tlj. We reproduce well the data 
in two different cases Xq = 100 and Xq = 0.032. The correct description of the px 
spectra in a broad range of rapidities is usually understood as an indication for a 
high level of thermalization of the entire fireball \TT\. After appropriate rescaling 
of the initial central energy density, the spectra are almost insensitive to initial 
values of anisotropy and to the length of the anisotropic stage (note a relatively 
large value of in Fig. 19. 3p . The small discrepancies in normalizations between 
different rapidity windows are mainly due to discrepancies between theoretical and 
experimental dNch/drjps profiles (compare with the part (c) of Fig. 19. ip . 
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Figure 9.3: Transverse momentum spectra of tt~^ (left panel) and (right panel) for 
the centrality class c = — 5% (5 = 2.23 fm) and for different rapidity windows. The 
model results are obtained with the time-scale parameter Teq = 1.0 fm. The results are 
obtained for two values of the initial anisotropy parameter: xq = 100 (dashed blue lines) 
and xo = 0.032 (dotted green lines). The model results are compared to the experimental 
data from BRAHMS [74J. The spectra for different rapidity windows are successively 
rescaled down by factor 0.1. 
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Figure 9.4: Left panels: Transverse momentum spectra of vr^ (top), (middle) and p 
(bottom) for centralities c = — 5% (6 = 2.26 fm) and c = 20 — 30% (6 = 7.16 fm) com- 
pared to the experimental data from PHENIX [102J. Right panels: Transverse- momentum 
spectra of A" (top), E~ (middle) and Q,~ + (bottom) for centralities c = — 5% and 
c = 20 — 40% (6 = 7.84 fm) compared to the experimental data from STAR |lU3j . The 
spectra for c = 20 — 30% and c = 20 — 40% are scaled down by factor 0.1. The results are 
obtained with time-scale parameter Tgq = 1.0 fm. Errors are statistical. The A feed-down 
corrections for proton spectra have be applied. The A spectra were feed-down corrected 
for S and Q. 
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In the left part of Fig. 19.41 we present the transverse-momentum spectra of posi- 
tive pions (top), kaons (middle), and protons (bottom) at midrapidity (|y| < 1) for 
the time-scale parameter Tgq = 1.0 fm and for centralities c = — 5% {b = 2.26 fm) 
and c = 20 — 30% (6 = 7.16 fm). The model results are compared to the data 
from PHENIX |102] . We can see that the model reproduces very well the slope of 
the spectra of ii'+'s and protons up to pt = 3 GeV. The model pion spectra for 
c = 20 — 30% are underestimated for pt > 2.5 GeV. This is an expected feature, 
since hard particles require more rescatterings to thermalize and this region is out 
of applicability of hydrodynamics. 

The correct normalization of spectra for different hadron species is assured by 
introducing thermodynamic parameters at freeze-out, which were obtained in the 
framework of a thermal model using a complete treatment of resonances [91]. The 
reproduction of the convex curvature of the pion spectra at low px is a consequence 
of feeding from resonance decays which are implemented in THEMINATOR [Ml ES]. 

On the right- hand- side of Fig. 19.41 we present the transverse-momentum spectra 
of hyperons for centralities c = — 5% (6 = 2.26 fm) and c = 20 — 40% {b = 7.84 fm) 
compared to the data from STAR [103] . Although our freeze-out temperature is a 
bit lower than that used in thermal models, the slopes of the hyperon spectra are 
reproduced very well. The normalization of A's is too small, the effect that may 
indicate their higher freeze-out temperature. Fig. 19.41 shows again a small sensitivity 
of the model predictions to initial anisotropy, which is consistent with conclusions 
based on Fig. 19. 3[ 

One would naively expect that the increase (decrease) of the pressure ratio P±/P\\ 
in the case xq = 100 (xq = 0.032), compared to the perfect-fiuid case, should 
lead to slower (faster) cooling of the system in the transverse direction. However 
the renormalization of the initial energy density (according to Table 19.11) in order 
to obtain the same final multiplicity leads to a reduction (amplification) of the 
transverse flow produced during expansion in the anisotropic phase [113j . This 
explicitly confirms the phenomenon of universality of the flow, predicting that the 
overall growth of the flow is the same regardless of the pressure anisotropy of the 
energy- momentum tensor ^106j . 



9.1.2 Anisotropic flows 

The significant amount of anisotropic flows is built at the very early stages of the 
evolution of matter, when the gradients of pressures are the largest due to the large 
spacial anisotropy of the source. The inclusion of early non-equilibrium stages of 
the collisions may significantly change the momentum anisotropy observed in the 
final spectra, making coefficients sensitive probes for the possible existence of 
non-equilibrium stages. 
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Figure 9.5: Transverse- momentum dependence of the elliptic flow coefficient of 7r'^+i^+ 
(left part) and protons (right part) calculated for the centrality c = 20 — 40% (6 = 7.84 fm) 
at midrapidity for the time-scale parameter Teq = 1.0 fm, and for two values of the initial 
anisotropy parameter. The results are compared to the PHENIX Collaboration data (red 
dots) [114] . The presented errors are statistical. Horizontal bars denote the p-r bins. 



Elliptic flow 

In the midrapidity region of non-central Au+Au collisions, the main information 
about the momentum anisotropy of the produced particles is contained in the 
coefficient. In the left part of Fig. 19.51 we show the elliptic flow of tt^ + as 
a function of the transverse momentum p-r- The model results obtained for the 
centrality c = 20—40 % (6 = 7.84 fm) and with the time-scale parameter Teq = 1.0 are 
compared to the data from PHENIX |114j . We observe that the model reproduces 
well the data up to = 1 GeV. The saturation for higher 'px is not reproduced, 
since it requires inclusion of shear viscosity effects [115] . Again, the model results 
weakly depend on the initial value of anisotropy, since the increase (decrease) of 
transverse pressure in the initial anisotropic stage is compensated by the decrease 
(increase) of the initial energy density. For completeness, in the right part of Fig. 19.51 
we show the p^-dependence of the elliptic flow of protons. The mass-splitting of V2 
is not reproduced. However, the elliptic flow of protons may be better described by 
introducing bulk viscosity }116] . 

In Fig. 19.61 we present the pr-integrated elliptic flow as a function of pseudo- 
rapidity. The model calculations are done for the centrality class c = 15 — 25%, 
three time-scale parameters, and three values of the initial anisotropy parameter. 
The model results are compared to the PHOBOS data [109] . We observe that the 
model overshoots the data by about 20% in the case of = 1 (solid black line). 
The peaked shape visible in the data is not well reproduced. 

We observe that the elliptic flow increases when the transverse pressure domi- 
nates. In this case, the elliptic flow is also steeper in pseudorapidity. It is easy to 
understand this behavior, since the reduction of longitudinal pressure leads to slower 
expansion of the fireball in the longitudinal direction. In the case of dominating lon- 
gitudinal pressure the behavior is opposite. The strength of this effect depends on 
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Figure 9.6: Pseudorapidity dependence of the pT-integrated elliptic flow of charged par- 
ticles calculated for centrality c = 15 — 25% for three time-scale parameters: Teq = 0.25 fm 
(a), Teq = 0.5 fm (b), Teq = 1.0 fm (c) and three values of initial anisotropy parameter. 
The results are compared to the PHOBOS Collaboration data (red dots) ^109j . The error 
bars show statistical errors. 



the length of the anisotropic stage. Thus, we conclude that the pr-integrated f 2 is 
sensitive to the early anisotropic stages. 

Directed flow 

When one considers non-central collisions of heavy-nuclei, one may observe non- 
zero directed flow in the region where rjps 7^ 0. This observable is sensitive to both 
transverse and longitudinal pressures, which makes it a perfect probe for measuring 
the early local anisotropics in momentum [113] . In Fig. 19. 71 we show the directed flow 
of charged particles for the two centrality bins c = — 5% (left part) and c = 5 — 40% 
{b = 6.79 fm, right part), and for the two time-scale parameters Teq = 0.25 fm (a) 
and Teq = 1.0 fm (b). The results are plotted together with the experimental points 
from STAR ^ and PHOBOS [11 Ij. 
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Figure 9.7: The pseudorapidity dependence of the directed flow of charged particles for 
two centrality bins c = — 5% (left part) and c = 5 — 40% (6 = 6.79 fm)(right part) for 
time-scale parameters Teq = 0.25 fm ((a) and (b)) and Teq = 1.0 fm ((c) and (d)) and 
three values of initial anisotropy parameter. The results are compared to experimental 
data from STAR (red dots) [llOj and PHOBOS (green squares) [lllj . Errors are statistical 
only. 



For different choices of our model parameters, the results describing vi are rather 
stable and consistent with the data. This leads to the conclusion that all the 
anisotropic stages described by ADHYDRO in this work are consistent with the Vi 
data. Moreover, the agreement with the data suggests the validity of the idea of a 
tilted initial source. 



9.1.3 HBT correlations 

The identical particle interferometry is a useful tool used in the analysis of the 
system sizes at freeze-out. In Fig. 19.81 we present the model results together with 
the STAR data |lU4j for the HBT radii: i?out (a), -Rside (b), and -Riong (c), all shown 
as functions of the pair total transverse momentum k^- Differences between the 
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Figure 9.8: The HBT correlation radii: i?out (a), Rside (b), i?iong (c) and ratio -Rout/-Rside 
(d) of positive pions as a function of total transverse momentum of the pair for centrality 
c = — 5% for time-scale parameter r^q = 1.0 fm and three values of initial anisotropy 
parameter. The results are compared to experimental data from STAR Collaboration (red 
dots) |lU4j . Errors bars reflect statistical errors. 



data and the model results are smaller than 10%. We observe that the slopes of the 
radii are quite well reproduced by our model. At large kr, the side radius, -Rside, 
is slightly smaller than that observed in the experiment, and the out radius, -Rout, 
is slightly larger. This is reflected in the overprediction of the -Rout /-Rside ratio by 
about 20%. The kx slope of -Rout /-Rside is not reproduced. 

For the all considered values of Xq, we obtain almost identical results. This is 
expected since the correct renormalization of the initial energy density results in 
the similar observed multiplicity and the flow on the freeze-out hypersurface (see 
our discussion in Sections 19.1.21 and 19. 2p . The similarities in the description of the 
correlation radii for all xq's mean that the shapes of the freeze-out hypersurfaces are 
also very much similar. The last point was already explicitly highlighted in Fig. 18.31 
in Chapter [HI 



80 



CHAPTER 9. NON-BOOST-INVARIANT DESCRIPTION 



9.2 Space dependent initial anisotropy 

In Ref. |117j it has been suggested that the system may be less equihbrated in 
forward and backward rapidities (i.e., in the fragmentation regions) than in the 
midrapidity region. On the other hand, in Ref. |118] it has been argued that the 
thermahzation may depend on the density of participants, since the degree of ther- 
mahzation depends on the number of colhsions between nucleons. Thus, one may 
expect that, in general, xq = xo{p{ri,yi±)), where xq ^ 1 at the edges of the system. 

Here we skip the possible dependence on centrality and check two scenarios: i) 
a constant initial anisotropy profile xo = 100 with Tgq = 1.0 fm, Si = 41.8 GeV/fm^, 
At] = 1.5, and cTj^ = 1.3, see Fig. 19.21 in Chapter 19.21 and ii) a space dependent 
initial anisotropy profile Xq = Xo{ri, yi±) where 

a;o(^,x^) = — , (9.2) 

with Teq = 1.0 fm, Ei = 73.8 GeV/fm^ At] = 1, and = 1.3. 

Fig. 19.91 shows the pseudorapidity dependence of charged particle multiplicity 
for the case i) (blue dashed line) and ii) (black dotted line). We observe that 
the distributions are similar despite different initial anisotropy profiles (note that 
in the case ii) At] is much smaller, since in the central rapidity region the matter 
evolves faster in the longitudinal direction than in the case i)). Moreover, moving 
to forward rapidities the anisotropy grows due to the ansatz (19. 2p . Thus, the tails of 
the pseudorapidity distribution are described more precisely. The results concerning 
the transverse momentum spectra are very similar in both i) and ii) cases. On the 
basis of conclusions from Section [9.1.11 this is an expected result. 




Figure 9.9: The pseudorapidity distributions of charged particles for the case i) (blue 
dashed line) and ii) (black dotted line). 

Interesting observations may be done for the f 2 coefficient at midrapidity. In the 
case i) the imposed initial anisotropy is large and constant in entire fireball. The 
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rescaling of the initial energy density cancels the effect of larger initial pressure in 
the transverse direction. Moreover the spatial anisotropy of the initially produced 
entropy in the case i), S(xo, cr(ro, x_l)), is the same as the spatial anisotropy of the 
initial entropy density profile. Hence the spatial anisotropy of the fireball transforms 
into the same momentum anisotropy as in the perfect-fluid case and we obtain similar 
results for V2 (see also Section [9.1.2p . In the case ii) the initial anisotropy is space 
dependent. In the center of the system it is almost isotropic due to high density of 
wounded nucleons. The anisotropy is much larger in the regions of low density of 
wounded nucleons, thus the spatial anisotropy of the initially produced entropy in 
the case ii) is lower than that in the case i). It results in a slightly lower spatial 
anisotropy of the whole fireball (see left part of Fig. 19.101) and in consequence it 
suppresses the generation of V2 (see right part of Fig. I9.10p . The effect of lowering of 
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Figure 9.10: Left panel: Time dependence of the spatial anisotropy of the fireball for 17 = 
for the case i) (blue dashed line) and ii) (black dotted line) with the entropy density as 
a weight function. Right panel: Transverse-momentum dependence of the elliptic flow 
coefficient of pions and kaons for the case i) (blue dashed line) and ii) (black dotted 
line). 

the f2 is clearly seen in the Fig. l9.11l where we present the pseudorapidity dependence 
of the pr-integrated V2- We observe better agreement of the calculations with the 
experimental data in the central rapidities in the case ii). However the model still 
does not reproduce the peaked shape seen in the data. The spatial dependence of 
the initial anisotropy has a small effect on the HBT radii. In Fig. 19.121 we present 
the HBT radii calculated for the case i) (blue dashed line) and ii) (black dotted 
line) . We observe a small improvement of -Rside which results in better agreement of 
-RoutZ-Rside in the case ii) as compared to the case i). 
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Figure 9.11: Pseudorapidity dependence of the py-integrated elliptic flow of charged 
particles calculated for the case i) (blue dashed line) and ii) (black dotted line). 
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Figure 9.12: The HBT correlation radii: i?out (a), -Rside (b), -Riong (c) and ratio 
-RoutZ-Rsidc (d) of positive pions as a function of total transverse momentum of the pair 
for centrality c = — 5% for the case i) (blue dashed line) and ii) (black dotted line). 



Chapter 10 

Summary and conclusions 



In this Thesis, a recently developed framework of highly-anisotropic and strongly- 
dissipative hydrodynamics - ADHYDRO - has been introduced and used to analyze 
the space-time evolution of matter produced in ultra-relativistic heavy-ion collisions. 
The main goal of this analysis was to study the effect of initial highly-anisotropic 
stages on the final soft hadronic obscrvables typically measured in the experiment. 
The study was done in the context of the heavy-ion measurements performed at 
RHIC (Relativistic Heavy Ion Collider) w Brookhaven National Laboratory. 

Starting from the general assumption about the form of the phase-space distribu- 
tion function, thermodynamic properties of locally anisotropic systems of particles 
have been studied and the form of the generalized equation of state has been for- 
mulated. The dynamic equations determining the evolution of a highly-anisotropic 
fluid have been introduced. The form of the entropy source related to the mech- 
anisms leading to thermalization of the system has been defined. In the simplest 
case of purely-longitudinal and boost-invariant expansion, different features of the 
model have been analyzed. In particular, for small deviations from equilibrium, the 
consistency of the ADHYDRO model with the Israel-Stewart theory has been demon- 
strated. 

Using the ADHYDRO model in the general (3+l)D and boost-invariant (2+l)D 
versions, the three possible scenarios of early stages of heavy-ion coUisions have 
been analyzed: i) the system initially oblate in momentum space, ii) the system 
initially prolate in momentum space, iii) the initially isotropic system. The results 
of the ADHYDRO model have been compared to the results obtained from the refer- 
ence, perfect-fiuid hydrodynamic model - LHYQUID. The results of the hydrodynamic 
models were coupled to the statistical Monte-Carlo model THERMINATOR which has 
allowed us to perform systematic study of soft-hadronic obscrvables: transverse- 
momentum spectra, directed and elliptic fiows, and the HBT radii. 

The following final conclusions may be drawn: 

a) All studied observables are almost insensitive to the initial anisotropic stage 
provided the initial conditions of the evolution are properly readjusted. 
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b) Complete thermalization of matter may take place only at the times of about 
1 fm/c. In this way the early thermalization puzzle may be circumvented. 

Finally, we have studied the systems where the initial anisotropy is spatially 
dependent. Our results suggest that the system may be highly anisotropic at the 
edges and simultaneously well thermalizcd in the center. It should be stressed, 
however, that the correct description of the data has been achieved only in the 
situations where the complete thermalization has been reached at about 1 fm/c. 



Appendix A 
Notation 



Here we gather the symbols used in the Thesis. 



Locally anisotropic systems of particles 



A_L,A|| transverse and longitudinal temperature, 

Pj_,P|| transverse and longitudinal pressure, 

m mass, 

mT,m± transverse mass, 

Pt,P± transverse momentum, 

Pl,P^ longitudinal momentum, 

X = (A^/A||) = ^ + 1 pressure anisotropy parameter, 

= ^ {l,Vx,Vy,Vz) four- velocity, 

y* = 7^ (t;^, 0, 0, 1) four- vector defining longitudinal direction, 

7 = (1 — Lorentz factor, 

= (1 ~ 'Vz)~^ Lorentz factor, 

p^ particle four-momentum, 

T^'^ energy-momentum tensor, 

S'^ entropy current, 

particle number current, 

e,a,n energy density, entropy density, and particle density 
in the anisotropic system, 

Eoq, o"cq. "-cq energy density, entropy density, and particle density 
in equilibrium, 

T temperature in the isotropic system, 

e statistics identifier (+1 for FD, -1 for BE), 

go degeneracy factor related to internal quantum 

numbers of particles, 

R pressure relaxation function, 

r{x) = tt'^ R{x) / {3go) rescaled pressure relaxation function, 

Cg sound velocity, 

Tc critical temperature, 
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Hydrodynamics 



entropy source, 

X and y components of the transverse vector u_l = {ux,Uy), 
longitudinal fluid rapidity, 
longitudinal proper time, 

rapidity, 

space-time rapidity, 

pressure tensor, 

viscous shear tensor, 

projection tensor, 

viscous bulk pressure, 

33 component of the shear tensor, 

shear relaxation time, 

shear viscosity, 

volume expansion rate, 

time-scale parameter, 

inverse of the relaxation time, 

initial proper time for starting hydrodynamics, 



Initial conditions 



p 


density of sources, 


p 


normalized density of sources. 




initial energy density at the center of the system 




in most central collisions, 


Pw 


wounded-nucleon density. 


Pb 


binary-collisions density, 


Pws 


Woods-Saxon nuclear density profile, 


K 


mixing factor between wounded-nucleon. 




and binary-collisions densities. 


T 


thickness function. 


fin 


total inelastic NN cross-section, 


po,RA,a 


parameters for Woods-Saxon nuclear proflle. 


b 


impact vector. 


6 = |b| 


impact parameter. 


fiv) 


initial longitudinal profile. 




parameters of the initial longitudinal proflle. 


6 


step function, 


rriN 


nucleon mass. 


Yb 


rapidity of the participant nucleon. 


rim 


range of rapidity correlations. 



"& = arctanhwz 
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Freeze-out 



9 

B,l3,S,C 
0p = aictaii{py / px) 

C,0,<t> 
A 

d(C,</.,0) 



freeze-out temperature, 

element of the three-dimensional freeze-out hypersurface, 
spin degeneracy factor, 
chemical potential, 

baryon, isospin, strange, and charm chemical potentials, 
baryon number, third component of isospin, strangeness, 

and charm, 

azimuthal angle of the transverse momentum, 
angles parametrizing freeze-out hypersurface, 

scale, 

distance from the point (r = tq, x = 0, y = 0, 77 = 0) 
to the hypersurface point with coordinates (C, (j), 6). 
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